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Abstract 

A  generalized  computer  model  for  determining  where 
to  target  conventional  weapons  against  a  non- reinforced 
concrete  runway  is  developed.  This  development  consists  of 
taking  the  overall  problem,  runway  destruction,  and  dividing 
it  into  a  number  of  identical,  stochastically  independent 
subproblems.  Then  one  of  these  subproble-ms  is  used  to 
determine  an  aim  point  strategy  which  reduces  the  number  of 
weapons  required  to  achieve  the  level  of  runway  destruction 
desired. 

In  order  to  solve  the  subproblem,  a  search  algorithm 
is  developed.  Along  with  this  search  algorithm,  two  methods 
for  approximating  the  level  of  destruction,  given  an  aim 
point  strategy,  are  investigated.  Finally,  the  entire  model 
is  verified  and  "validated"  by  a  number  of  statistical  tests 
Also  a  comparison  of  this  model  to  the  Airfield  Assessment 
Program,  an  Air  Force  Armament  Laboratory  model,  is  provided. 
Using  the  model’s  results  for  the  data  investigated,  an  equa 
tion  relating  circular  error  probable  and  minimum  launch 
window's  width  to  the  number  of  weapons  employed  is  derived. 
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A  GENERALIZED  COMPUTER  MODEL  FOR  THE 


TARGETING  OF  CONVENTIONAL  WEAPONS 
TO  DESTROY  A  RUNWAY 

I.  Introduction 


Background 

The  airfield  complex,  normally  attacked  in  hopes  of 
reducing  an  enemy's  air  capability,  has  been  a  popular 
target.  Of  primary  importance  on  the  complex  have  been 
the  aircraft  themselves  because  of  their  delicate 
structure,  fuel  content,  and  high- explosive  weapons. 

Two  factors  have  made  aircraft  a  less  opportune  target 
and  centered  the  attack  on  airfield  pavement  surfaces. 
First,  hardened  shelters  have  made  the  aircraft  harder 
to  destroy.  The  attacker  must  breach  or  substantially 
damage  the  shelter  to  reach  the  aircraft.  This  requires 
special  weapons  whose  accuracy  demands  aerial  maneuvers 
which  make  the  attacking  planes  vulnerable  to  ground 
defenses.  The  second  factor  has  been  the  increase  in 
aircraft  performance  and  the  consequent  intensity  with 
which  the  first  stages  of  a  war  will  be  fought. 

Attacking  planes  can  strike,  return  to  base,  rearm, 
and  be  over  the  target  again  in  an  extremely  short  time, 
often  shorter  than  that  required  to  repair  a  runway  and 
launch  interceptors.  Opposition  fighters  can  thus  be 
effectively  eliminated  from  the  first  stages  of  battle 
without  being  physically  damaged  by  denying  them  a  launch 
and  recovery  surface.  The  concept  that  the  runways.  .  . 
are  poor  targets  is  no  longer  valid  in  the  light  of  the 
pace  of  modern  warfare  and  the  capabilities  of  aircraft 
[Ref  7:  260] . 

Currently,  there  is  extensive  research  being  conducted 
in  developing  conventional  weapons  for  runway  destruction, 
but  there  seems  to  be  no  readily  available  targeting  scheme 
on  how  best  to  employ  these  weapons.  There  are  a  number  of 
computer  models  used  by  the  Air  Force  to  simulate  runway 
destruction;  for  example,  the  Airfield  Assessment  Program, 
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AAP  (Ref  3) ,  but  these  models  require  the  user  to  supply  the 
aim  points  he  desired.  The  models  then  determine,  through 
Monte  Carlo  simulation,  whether  the  runway  has  been  sufficiently 
destroyed  so  a  specified  minimum  launch  window  no  longer  exists 
(minimum  launch  window  being  the  minimum  length  and  width  of 
runway  required  for  takeoff  or  landing) .  The  main  drawback 
to  these  investigations  is  the  locations  of  the  aim  points 
and  the  number  of  weapons  targeted  per  aim  point  (referred  to 
for  the  rest  of  this  thesis  as  the  "aim  point  strategy")  are 
mostly  derived  from  the  experience  of  the  user  and  his  best 
guess  as  to  what  this  strategy  should  be.  If  one  could  deter¬ 
mine  an  aim  point  strategy  which  reduces  the  total  number  of 
weapons  employed  and  still  maintain  the  desired  level  of  run¬ 
way  destruction,  then  one  could  lower  the  overall  costs  with¬ 
out  lowering  the  effectiveness  (assuming  that  all  other  costs 
were  equal) . 

Problem 

The  problem  is  to  develop  a  computer  model  which  will 
determine  an  aim  point  strategy  that  achieves  the  desired 
level  of  runway  destruction  with  as  few  total  weapons  as 
possible. 

Assumptions 

1.  The  circular  error  probable,  CEP,  of  the  weapons 
employed  can  be  logically  described  by  a  circular  normal  proba¬ 
bility  distribution. 

2.  The  survivability  of  all  weapons  is  one.  As  long 
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as  the  weapon's  survivability  is  the  same  for  all  aim  points, 
this  assumption  will  only  reduce  the  total  number  of  weapons 
employed  and  will  not  affect  the  number  and  the  locations  of 
the  aim  points. 

3.  The  reliability  of  all  weapons  is  one. 

4.  Each  weapon  is  independently  targeted  so  that  its 
impact  point  does  not  depend  on  any  other  weapon's  impact 
point  or  aim  point. 

5.  Only  one  type  of  weapon  can  be  used  to  destroy  the 
runway.  Thus,  the  CEP  ar  l  the  effective  explosive  weight  will 
remain  constant  throughout  an  investigation  (investigation 
being  the  determination  of  a  good  aim  point  strategy  for  a 
specific  set  of  inputs) . 

6.  The  time  frame  between  the  launching  of  the  first 
weapon  and  the  launching  of  the  last  weapon  is  too  short  to 
update  aim  points  through  damage  assessment. 

7.  Any  minimum  launch  window  must  be  totally  composed 
of  actual  runway. 

8.  If  a  minimum  launch  window  exists,  it  is  assumed 
to  be  accessible  to  an  aircraft,  even  though  this  may  not  be 
true  (Ref  3) . 

9.  The  runway  destruction  criteria  is  a  given  confi¬ 
dence  level  that  no  minimum  launch  window  will  remain  undam¬ 
aged  after  the  attack. 

Scope 

Aside  from  the  limitations  inherent  in  the  assumptions, 
the  model  is  designed  for  investigations  using  conventional 
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weapons  in  destroying  non-reinforced  concrete  runways.1  Yet 
the  model  is  general  enough  to  investigate  any  size  non- 
reinforced  concrete  runway  and  any  size  minimum  launch  window. 

The  model  was  used  to  investigate  destroying  a  runway 
8000  feet  long  and  150  feet  wide  (the  dimensions  considered 
typical  of  a  North  Atlantic  Treaty  Organization  runway  (Ref 
3))  and  a  minimum  launch  window  2000  feet  in  length  fa  length 
which  is  slightly  shorter  than  a  typical  fighter's  takeoff 
roll  (Ref  8)).  This  minimum  launch  window's  length  was 
coupled  with  a  width  of  either  50  feet  or  100  feet.  With 
these  sets  of  dimensions,  effective  explosive  weights  of  250, 
500,  and  1000  pounds  and  CEPs  of  20,  50,  and  100  feet  were 
investigated.  The  runway  destruction  criterion  for  ail  invest- 
gations  was  maintained  at  a  level  of  .80. 

Approach  and  Presentation 

The  approach  used  in  this  thesis  is  first  to  reduce 
the  two  dimensional  problem  to  a  set  of  independent,  identical 
one  dimensional  subproblems  (Chapter  II).  Second,  for  a  given 
aim  point  strategy,  two  methods  for  approximating  the  destruc¬ 
tion  level  are  developed  (Chapter  III).  Third,  a  search  algor¬ 
ithm  which  determines  a  good  aim  point  strategy  is  developed 
(Chapter  IV).  Next,  the  two  approximation  methods  and  the 
algorithm  are  combined  for  verification  and  validation  (Chap¬ 
ter  V).  Finally,  a  discussion  of  the  conclusions  and  recom¬ 
mendations  for  further  investigation  is  presented  (Chapter  VI). 

1Atypical  of  current  runways  (see  Recommendation  5) 
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II.  Reduction  of  the  Two  Dimensional 


Problem  to  a  One  Dimensional  Problem 
Introduction 

As  was  pointed  out  in  Chapter  I ,  the  approach  used  in 
this  thesis  is  to  reduce  a  two  dimensional  problem  (the  two 
dimensions  being  length  and  width)  to  a  one  dimensional  prob¬ 
lem  (just  width).  In  order  to  accomplish  this  reduction,  the 
independent  variables  (the  "givens"  for  the  problem)  must 
first  be  transformed  into  a  new  set  of  variables.  Next,  the 
overall  runway  problem  will  be  subdivided  into  a  number  of 
like,  smaller,  stochastically  independent  subproblems  (the 
runway  "cuts”).  Since  the  subproblems  are  alike  and  stochas¬ 
tically  independent,  by  minimizing  the  number  of  weapons  re¬ 
quired  to  accomplish  any  one  cut,  the  minimum  number  of 
weapons  for  any  cut  is  then  determined.  Finally,  the  two 
dimensional  subproblem  (any  one  cut)  is  transformed  into  a 
one  dimensional  subproblem. 

Transformation  of  the 
Independent  Variable's 

The  five  independent  variables  the  user  is  required  to 
provide  are: 

1.  The  weapon's  effective  explosive  weight,  W,  in 
pounds  TNT 

2.  The  weapon's  circular  error  probable,  CEP,  in  feet 
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3.  The  actual  runway  length,  X,  and  width,  Y,  in  feet 

4.  The  length,  U,  and  width,  V,  of  the  minimum  launch 


window 

5.  The  runway  destruction  criterion,  P  (a  real  number 
between  zero  and  one  which  represents  the  desired 
confidence  level  that  no  minimum  launch  window 
remains  undamaged  after  the  attack) 

Four  of  these  independent  variables  (W,  CEP,  X  and  Y,  and  P) 
are  converted  into  a  new  set  of  independent  variables  (damage 
radius,  standard  deviation,  effective  runway  length  and  effec¬ 
tive  runway  width,  and  the  cut  destruction  criterion). 

Effective  Explosive  Weight.  The  effective  explosive 
weight  is  converted  into  the  damage  radius,  R,  this  explosive 
would  cause  in  non-reinforced  concrete  runways.  This  conver¬ 
sion  is  accomplished  by  determining  the  optimal  depth  of 
burst  (DOB)  for  a  specific  weapon,  the  damage  radius  caused 
by  this  specific  weapon  at  its  optimal  DOB,  and  a  scaling 
factor.  Coupling  these  three  items,  one  can  derive  the  fol¬ 
lowing  equation  (see  Appendix  A  for  the  derivation) : 

R  =  (3.54)(W)1/3  (2-1) 

Circular  Error  Probable.  The  second  independent  vari¬ 
able  transformed  is  the  circular  error  probable,  CEP,  for  the 
weapon.  CEP  is  the  radius  of  a  circle  which  encloses  50  per¬ 
cent  of  the  weapons'  impact  points  given  that  the  weapons  are 
targeted  for  the  center  of  this  circle  and  there  are  no  biasing 
errors.  For  example,  if  ten  weapons  with  identical  CEPs  of 
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50  feet  are  targeted  at  an  aim  point,  on  the  average  five  of 
these  weapons  will  impact  within  50  feet  of  that  aim  point. 

The  approach  used  by  this  thesis  required  that  CEP  be  converted 
into  the  standard  deviation  (o)  of  the  weapon.  The  conversion 
factor  is 


CEP 

(2  In  2)1/2 


(2-2) 


(see  Appendix  B  for  the  derivation) . 

Actual  Runway  Length  and  Width.  The  third  transforma¬ 
tion  is  from  actual  runway  length  and  width  to  effective  run¬ 
way  length  and  width.  The  effective  runway  dimensions  are 
defined  to  be  those  dimensions  within  which  a  weapon  can  impact 
and  cause  damage  to  the  actual  runway.  The  effective  runway’s 
dimensions  will  be  greater  than  those  of  the  actual  runway 
because  a  weapon  can  impact  off  the  actual  runway  up  to  a  dis¬ 
tance  of  R  and  still  cause  damage  to  the  runway.  The  effective 
runway  length,  EX,  is 


EX  =  X  +  2R  (2-3) 

and  the  effective  runway  width,  EY,  is 

EY  =  Y  +  2R  (2-4) 

Using  these  dimensions,  EX  by  EY,  there  are  four  areas 
in  which  a  weapon  could  impact  and  still  not  cause  damage  to 
the  actual  runway.  These  four  areas  are  near  the  corners  of 
the  rectangle  formed  by  EX  and  EY.  Even  though  these  areas 
do  exist,  when  compared  to  the  entire  area  under  investigation, 
they  are  insignificant  and  will  be  ignored. 
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Runway  Destruction  Criterion.  The  final  transformation 
is  to  convert  the  runway  destruction  criterion,  P,  to  the  des¬ 
truction  criterion,  PC,  for  one  of  the  subproblems.  As  a 
first  cut,  it  is  assumed  the  subproblems  are  identical  and 
stochastically  independent,  thus 

PC  =  (P)1/N  (2-5) 

where  N  is  the  number  of  subproblems.  (Note:  when  one  uses 
whole  numbers  of  weapons  to  achieve  PC  and  requires  each  sub¬ 
problem  to  be  identical,  there  may  be  some  "overkill"  when 
the  subproblems  are  combined  into  the  final  solution.) 

Subdividing  the  Overall 
Runway  Problem 

One  now  has  the  following  set  of  independent  variables: 

1.  The  damage  radius,  R,  of  the  weapon 

2.  The  standard  deviation,  a,  of  the  weapon  impact 
pattern 

3.  The  effective  runway  length,  EX,  and  the  effective 
runway  width,  EY 

4.  The  length,  U,  and  width,  V,  of  the  minimum  launch 
window 

5.  The  destruction  criterion  for  each  of  the  subprob¬ 
lems  ,  PC 

The  approach  used  in  this  thesis  is  to  cut  the  runway 
across  the  effective  width  in  enough  locations  so  lengthwise 
the  runway  destruction  criterion  is  achieved.  The  size  of  each 
of  these  weapon  impact  areas,  cut,  will  be  6a  by  EY  (see  Fig  2). 

6a  is  chosen  as  the  cut’s  width  to  help  simplify  the 
subproblem.  If  a  weapon  is  targeted  for  the  centerline  of  a 

cut,  one  is  99  percent  confident  that  the  weapon  will  impact 
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within  3 a  of  the  centerline.  Thus,  the  probability  that  this 
weapon  impacts  further  than  3a  from  the  centerline  is  extremely 
small,  and  will  be  assumed  to  be  zero.  In  addition  to  the  cut 
dimensions,  one  must  determine  the  minimum  takeoff  distance 
along  the  runway's  length.  If  an  aircraft  takes  off  across 
the  runway  rather  than  along  the  axis  of  the  runway's  length, 
it  would  require  less  of  the  actual  runway's  length.  By 
knowing  U  and  Y,  the  minimum  takeoff  distance  (along  the  length 
axis  of  the  runway) ,  D,  is  (see  Fig  1)  : 

D  =  (U2  -  Y2)1/2  (2-6) 

Even  though  this  type  of  takeoff  violates  Assumption  7  (Chap¬ 
ter  I),  D  is  less  than  any  minimum  takeoff  distance  required 
under  Assumption  7.  Thus,  if  the  runway  destruction  criteria 
is  met  using  D,  it  is  also  met  under  Assumption  7. 

D  can  now  be  used  to  determine  the  interval  between 
any  two  adjacent  cuts  on  the  runway.  Coupling  the  dimensions 
of  the  cut  with  D,  the  interval  between  the  centerlines  of  any 
two  adjacent  cuts,  B,  is  (see  Fig  2): 

B  =  D  -  6a  (2-7) 

(An  implicit  assumption  in  this  approach  is  that  no  two  cuts 
will  overlap.  If  they  do  overlap,  the  cuts  are  no  longer 
stochastically  independent,  and  PC  must  be  computed  using 
conditional  probabilities.)  B  is  appropriate  for  any  two 
adjacent  cuts,  but  it  is  inappropriate  for  the  leading  edge 
of  the  actual  runway  and  the  trailing  edge  of  the  first  cut. 


Leading  will  refer  to  the  left  extreme  of  the  object  associ¬ 
ated  with  it;  for  instance,  the  actual  runway.  Trailing  will 
indicate  the  right  extreme  of  the  object  it  is  associated  with; 
in  this  case,  the  first  cut  (as  indicated  in  Fig  3).  The 
interval  between  the  leading  edge  of  the  actual  runway  and 
the  centerline  of  the  first  cut  (the  first  interval)  is  equal 
to  D  -  3o.  Thus,  the  total  number  of  cuts,  N,  to  be  made  in 
the  runway  is 

N  =  [x-  CD~.3g)  ]  +  l  (2-8) 

where  the  symbol  [a]  indicates  the  largest  integer  value  less 
than  or  equal  to  a.  N  must  be  an  integer  value  because  any 
remainder  can  be  viewed  as  if  it  is  the  distance  from  the 
leading  edge  of  the  last  cut  to  the  trailing  edge  of  the 
actual  runway.  This  distance  must  be  less  than  D,  and  thus 
it  is  less  than  the  length  of  the  minimum  launch  window.  The 
+1  on  the  right-hand  side  of  Eq  (2-8)  accounts  for  the  first 
cut  which  is  being  subtracted  from  X  in  order  to  determine 
the  number  of  remaining  cuts  (D  is  the  interval  between  any 
two  adjacent  cuts,  and  this  interval  is  less  than  the  first 
interval) . 

If  the  N  cuts  are  positioned  at  their  proper  intervals 
along  the  runway,  any  minimum  launch  window  must  pass  through 
at  least  one  of  these  cuts  (Assumption  7  forces  this  to  be 
true).  Also,  the  angular  difference  between  a  minimum  launch 
window  along  the  length  axis  of  the  runway  and  a  minimum  launch 
window  along  any  axis  which  satisfied  Assumption  7  (an  angle 
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Leading  Edge 
of  the 

Actual  Runway 


I 


1 

1 

<  30 

1 

1 

1 

D 


Trailing  Edge 
of  the 
First  Cut 


Fig  3.  Interval  from  the  Leading  Edge 
of  the  Actual  Runway  to  the 
Trailing  Edge  of  the  First  Cut 


less  than  angle  b  in  Fig  1)  is  very  small.  One  can  thus 
assume  that  if  PC  is  met  for  an  axis  along  the  length  of  the 
runway,  PC  is  effectively  met  for  any  axis.  Furthermore,  if 
all  cuts  are  achieved  at  a  level  of  PC,  the  actual  runway 
destruction  criteria  is  met;  this  discussion  should  justify 
Eq  (2-5). 

Converting  the  Subproblem 
to  One  Dimension 

By  combining  two  previously  discussed  ideas,  the  two 
dimensional  subproblem  can  be  transformed  into  a  one  dimen¬ 
sional  subproblem.  These  two  ideas  are: 

1.  If  PC  is  met  for  any  axis  along  the  length  of  the 
runway,  then  PC  is  effectively  met  for  any  axis. 

2.  Any  minimum  launch  window  must  pass  entirely 
through  one  cut. 

Thus,  all  impact  points  within  a  cut  may  be  viewed  as  if  they 
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are  on  the  cut's  centerline.  Also,  due  to  the  assumption  used 
in  choosing  the  cut's  width,  the  circular  normal  probability 
distribution  can  be  viewed  as  a  normal  probability  distribu¬ 
tion  in  one  dimension.  Finally,  since  all  cuts  are  identical 
in  size  and  PC,  one  only  needs  to  determine  the  optimal  aim 
point  strategy  for  a  single  cut  to  have  the  optimal  aim  point 
strategy  for  all  cuts  (the  only  real  difference  is  in  the 
locations  of  the  cuts'  centerline).  Furthermore,  due  to  the 
exponential  form  of  the  cumulative  probability  distribution 
for  any  PC,  fewer  weapons  are  required  to  satisfy  P  if  all 
cuts  have  the  same  probability  of  destruction.  Thus,  the 
minimum  total  number  of  weapons  required  to  achieve  P  is  N 
times  the  number  of  weapons  needed  to  make  any  one  cut  at  a 
confidence  level  of  PC. 

Summary 

The  overall  two  dimensional  problem  has  been  converted 
to  a  one  dimensional  subproblem.  This  was  accomplished  by 
analytically  transforming  four  of  the  five  independent  vari¬ 
ables  into  new  independent  variables.  This  new  set  of  inde¬ 
pendent  variables  was  then  used  to  logically  divide  the  over¬ 
all  problem  into  a  number  of  like,  stochastically  independent 
subproblems  in  two  dimensions.  Finally,  due  to  the  stated 
assumptions  and  the  formulation  of  the  suhproblems,  each  sub¬ 
problem  can  be  viewed  in  one  dimension  and  by  solving  only  one 
of  the  subproblems,  the  entire  problem  is  solved. 
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III.  Methods  for  Approximating  the 


Probability  of  Achieving  a  Cut 

Introduction 

Using  the  set  of  model  parameters  and  the  simplifica¬ 
tion  of  the  overall  problem  discussed  in  Chapter  II,  the  one 
dimensional  subproblem  can  be  expressed  as  follows: 

Determine  the  aim  point  strategy  which 
Minimizes  the  number  of  weapons  required 
Subject  to  no  section,  wider  than  the  minimum  (3-1) 
launch  window’s  width,  remains  un¬ 
damaged  with  a  probability  of  at 
least  PC 

In  order  to  develop  an  algorithm  which  approximates 
the  solution  to  this  subproblem,  a  building  block  approach 
was  used.  The  first  phase  was  to  determine  an  adequate  approx¬ 
imation  for  the  probability  of  making  a  cut,  PC,  given  an  aim 
point  strategy.  This  required  choosing  the  approximation 
methods  and  then  adapting  the  subproblem  to  the  approximations 
used. 

Approximation  Methods  for  a 
Given  Aim  Point  Strategy 

Three  methods  were  investigated.  The  first  method  in¬ 
corporates  order  statistics.  The  second  uses  the  event- 
composition  method.  And  the  third  method  involves  a  Monte 


14 


Carlo  simulation. 


Order  Statistics  Method.  Under  this  method,  the  con¬ 
straint  can  be  formulated  as  follows: 

1.  Let  H  be  the  probability  distribution  of  the  weapons' 
impact  locations  on  the  effective  runway  for  a 
given  aim  point  strategy. 


Let 

h^  be  the 

realization  of 

the 

ith 

order  statis 

tic 

,  where  h^ 

is  one  edge  of 

EY; 

h(N) 

is  the 

other  edge  of  EY ;  and  h^  ,  i  =  1,2,...,N-1,  are 
the  impact  points  themselves. 

3.  Let  G  be  the  probability  distribution  of  the 
maximum  difference  between  h^  and  h^+^ 

Then 

P(g  <  V  +  2R)  >  PC  (3-2) 

In  order  to  use  this  method,  one  first  needs  to  approximate 
the  probability  distribution  of  H.  But  the  distribution  of 
H,  in  most  cases,  is  the  joint  distribution  of  more  than  one 
normal  probability  distribution.  Coupling  this  fact  with  the 
requirement  that  h^  could  be  from  any  of  the  normal  distri¬ 
butions  about  their  aim  points  for  the  given  strategy,  the 
order  statistic  method  was  considered  not  promising  as  a 
model  for  estimating  PC. 

Event-Composition  Method.  Under  the  event-composition 
method,  the  event  of  interest,  the  achieving  of  a  cut,  is 
expressed  "as  a  composition  (unions  and/or  intersections)  of 
two  or  more  other  events  [Ref  9:  42]."  This  method  requires 
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transforming  a  continuous  probability  distribution,  I  (the 
probability  distribution  of  achieving  a  cut  with  a  given  aim 
point  strategy),  into  a  set  of  discrete  probability  subsets. 

Monte  Carlo  Simulation.  With  the  Monte  Carlo  simulation, 
the  impact  points  of  the  weapons  are  generated  using  a  normal 
random  generator  with  a  mean  of  zero  and  a  standard  deviation 
of  one.  By  coupling  these  random  numbers  with  the  aim  point 
strategy  and  the  standard  deviation  of  the  weapons'  impact 
points,  one  can  produce  a  set  of  random  variates  which  will 
duplicate  an  expected  attack.  The  set  of  variates  can  then 
be  tested  to  see  if  this  realization  of  the  aim  point  strategy 
cut  the  runway. 

Description  of  the 
Event-Composition  Method 

Rather  than  approximating  PC  for  a  given  aim  point 
strategy,  the  complement  of  PC,  FU,  was  chosen.  In  order  to 
use  the  event- composition  method,  one  must  first  define  the 
events  which  comprise  the  overall  event  of  interest. 

1.  E  is  the  event  that  at  least  one  occurs  in  an 
attack  (by  definition  P(E)  =  FC)  . 

2.  is  the  event  that  no  weapon  damages  the  it}l 
minimum  launch  window's  width  (this  event  occurs 
when  no  weapon  impacts  in  the  it^1  grouping)  . 

3.  is  the  event  that  no  weapon  targeted  for  the 

j1*1  aim  point  impacts  in  the  k**1  section  (length  Z) 
th 

of  the  i  grouping.  (For  further  understanding 
of  these  events,  refer  to  Fig  4) 
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Given:  one  aim  point  (j  =  1) 

four  sections  per  group  (k  =  1,2, 3, 4) 
two  groupings  (i  =  1,2) 


where  the  small  a-,.,  1.,  e 

jki*  i* 

presentation  of  events 
respectively. 


represent  a  graphic 
L . ,  and  E 

l  ’ 


Fig  4.  A  Graphic  Representation  of 
Events  Ajki,  L^,  and  E 

The  Numbers  of  the  Events  Which  Comprise  E.  Prior  to 
determining  the  probability  of  a  specific  Ajki,  P(Ajki),  one 
must  first  know  the  length  Z,  part  of  the  definition  of  Ajki. 

In  determining  Z  for  a  given  aim  point  strategy,  there  is  an 
inverse  relationship  between  the  length  of  Z  and  the  accuracy 
of  this  method  in  approximating  the  actual  continuous  situation, 
i.e.  as  Z  becomes  longer,  this  method  becomes  less  accurate. 

Yet  as  Z  decreases,  user's  inputs  remaining  constant,  the  com¬ 
putation  time  increases. 

It  is  felt  the  user's  minimum  launch  window's  width,  V, 
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is  a  point  estimate  of  the  actual  required  minimum  launch 
window's  width.  This  being  true,  the  results  for  any  con¬ 
tinuous  case  based  on  the  user's  inputs  would,  in  itself,  be 
an  approximation.  As  a  first  decision  was  set  to  .IV. 

This  Z  value,  Z^,  is  used  for  all  sections  within  a  grouping 
except  for  the  last  section  within  each  grouping. 

In  order  to  cause  any  estimate  of  failing  to  cut  the 
runway  to  be  higher  than  is  true,  or  "optimistic",  the  follow 
ing  formulas  were  chosen  to  determine  the  number  of  sections 
per  width,  NA;  the  number  of  sections  per  grouping,  NL;  and 
the  number  of  groupings  per  width,  NW. 

NA  =  +  1]  (3-3) 

where  Y  is  the  actual  runway  width  and  R  is  the  damage  radius 

NL  =  [Y-*^—  +  1]  (3-4) 

where  F  is 

F  -  -  [£])R  (3-5) 

NW  =  NA  -  NL  +  1  (3-6) 

The  +1  on  the  right-hand  sides  of  Eqs  (3-3)  and  (3-4) 
allows  one  to  incorporate  the  fraction,  F,  into  each  grouping 
for  a  specific  L^.  The  last  section  in  each  grouping  (Z2  in 
size),  except  for  the  last  grouping,  is  actually  shorter  than 
the  other  sections  within  the  grouping.  The  last  section 
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the  last  grouping  (Z^  in  size)  not  only  incorporates  the 
fraction,  F^,  it  also  includes  the  remaining  effective  runway 
width  not  accounted  for  to  this  point.  Thus, 


Zl  £or  ^j,i  and  k  =  1,...,NL-1 
Z  for  V 

2  j;  i  =  l , . . . NW- 1 ;  k=NL 

Z 3  £or  Vj;  i=NW;  k=NL 


(3-7) 


With  the  numbers  of  events  which  comprise  the  subproblem,  one 
can  use  this  information  in  determining  the  probability  of  E, 
P(E).  The  approach  to  accomplishing  this  task  will  start  with 
Ajki>  use  the  A^^'s  to  build  the  I^'s,  and  then  use  the  L^'s 
to  build  E. 

Probability  of  Event  A^^.  Each  Aj ^  within  grouping  i 
are  designed  to  be  mutually  exclusive  (see  Fig  4).  Yet,  they 
are  not  independent  of  each  other.  If  we  know  that  A^^  occurs, 
then  any  revised  probability  of  A^^  occuring  will  decrease. 
Assuming  there  are  N  aim  points  (j  =  1,...,N)  and  one  weapon 
per  aim  point,  the  probability  of  A^^,  p(A^i)  is 


P<Ajki>  ’  1  -  jlj  PtAjlci’  Vk,i 


(3-8) 


where  Aj^  is  the  complement  of  Aj^. 

Probability  of  Event  L^.  By  definition  only  occurs 
when  all  A^^  within  grouping  i  occur.  This  can  be  expressed 


as  follows: 


Ali  A  A2i  A",A  ^SlL  ,  i 


(3-9) 
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or 


Li  A  li  U  A  2i  U*  *  *U  SNL,i 


(3-10) 


Since  all  within  a  grouping  i  are  mutually  exclusive,  the 
probability  of  L^,  P(L^)  is  (Ref  9:  44) 


NL 


P(L 


)  =  1  -  Z  P(I  k  •) 
1  k=l  1»x»1 


(3-11) 


(Note:  this  refers  to  one  weapon  at  aim  point  1) 

By  expanding  the  number  of  aim  points  to  N  and  realizing  that 
events  from  different  aim  points  are  independent, 

N  NL 

P (L - )  =  n  (1  -  Z  P(A.V .))  (3-12) 

1  j=l  k=l  JK1 

(Note:  still  only  one  weapon  per  aim  point) 

Probability  of  Event  E.  As  previously  defined,  E  is 

the  union  of  all  L^s,  thus 


E  =  L1  U  L 


The  probability  of  event  E, 


NW 

P(E)  =  Z  P(L  • ) 
i  =  l  1 


NW-  2  NW- 1  NW 
♦  Z  Z  Z  P(L 
a-1  b=2  c=3  a 

a<b<c 


2  U • • • U  LNW 


P (E)  is  (Ref  9:  46) 
NW  - 1  NW 

-  Z  Z  P(L  A  Lb) 
a=l  b=  2  a  D 
a<b 


A  Lb  A  Lc)  -  . . . +. 


-  (-1)NW  P(LX  A  L2  A. ..A  Lnw) 
Incorporating  P(L^),  Eq  (3-12): 


(3-13) 


+ 


(3-14) 
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NW  N  NL 

P(E)  =  Z  {  n  (1-  Z  P(K..  .))}  - 
i=l  j=l  k=l  JK1 

NW- 1  NW  N  NL  MIN(NL,b-a) 

z  z  {  n  (l  -  z  P(I.kh)  -  z  P(ft,kJ)) 

a=l  b=2  j=l  k=l  k=l  JKa 

a<b 

NW  -  2  NW- 1  NW  N  NL 

+  I  Z  Z  {  n  (1  -  Z  P(E..  )  - 

a=l  b=2  c=3  j=l  k=l  J KC 

a<b<c 

MIN(NL,c-b)  _  MIN(NL.b-a) 

Z  P  (A  •  vh)  -  Z  P(^.;k,,))}  + 

k=l  Jlcb  k=l  JKa 

MW  N  NL  . 

-  (-1)  {.n  (l  -  f/^.k.NW^  '  p^j,i»NW-i)  "  **•  ' 

J  (3-15) 

where  the  symbol  MIN(a,b)  means  the  minimum  of  a  and  b. 

The  only  feature  which  P(E)  of  Eq  (3-15)  does  not  incor¬ 
porate  is  the  number  of  weapons  targeted  for  a  specific  aim 
point  (other  than  the  trivial  case  of  one  weapon  per  aim  point) . 
By  letting  be  the  number  of  weapons  targeted  for  the  aim 
point,  P(L^)  becomes 


NL  T. 

P(L  )  =  (1  -  Z  PCS...))  J 
i  k=l  JKi 


(3-16) 


and  thus, 


NW  N  NL  T . 

P(E)  =  Z  {  n  (1  -  Z  P(S...))  J} 
i=l  j=l  k=l  JK1 


NW-1  NW  N  NL 

■  Z  Z  {  n  (1  -  Z  P(X,kh) 
a*l  b=2  j=l  k=l  JKD 

a<b 


MIN(NL,b-a)  T. 

jkb3  ‘  P(*jka)5  } 
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-  (~ij  ”  ln  (l  -  z  P(£.k  NW)  - 
j=l  k=l  JK,NW 


p(^j ,l,NW-l^ 


-  P(^  i  ,))  j> 

J  t  x  f  A 


(3-17) 


(See  Appendix  C,  Subroutine  UNION  for  the  coding  of  Eq  (3-17).) 

Solving  for  PfS^,.).  With  P(E)  as  formulated,  the  only 
J  K1 

probabilities  which  must  be  determined  to  solve  Eq  (3-17)  are 
the  P (^j lci^  * s  *  s^nce  distribution  of  the  weapons'  impact 
points  around  their  aim  points  is  a  normal  probability  distri¬ 
bution,  one  needs  to: 

1.  Translate  EY  such  that  the  location  of  aim  point  j 
is  assigned  a  value  of  zero. 

2.  Convert  the  translated  leading  edge,  d,  and  the 

translated  trailing  edge,  f,  of  each  into 

standard  deviations  about  aim  point  j.  In  a  majority 
of  cases,  the  trailing  edge  of  a^k^  is  the  leading 
edge  of  aj>kn>i. 

3.  Solve  the  following  integral: 


1  1/0  12 

P(d  <  x  <  f)  =  -  f  exp  -  i(x  jdx 

/27  d/a  z 


(3-18) 


where  P(d  <_  x  <_  f)  is  the  probability  a  weapon 
targeted  for  aim  point  j  will  impact  such  that  A... 

J  K1 

will  occur.  (See  Appendix  C,  Subroutine  PROCAL 
for  the  coding  of  this  transformation) 
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Description  of  the  Monte 
Carlo  Simulation  Method 

As  with  the  event-composition  method,  the  Monte  Carlo 
simulation  method  approximates  for  a  given  aim  point  stra¬ 
tegy.  As  was  explained  in  the  paragraph  Monte  Carlo  Simulation, 
the  set  of  random  variates  (an  "attack")  either  cuts  the  run¬ 
way  or  it  does  not  cut  the  runway.  If  one  simulates  M, 
identical  and  independent,  "attacks",  then  this  method  becomes 
a  binomial  experiment  (Ref  9:  72).  And 

as  a  rough  rule,  you  can  assume  that  the  distribution 

A  A 

of  p  (p,  the  estimate  of  PTT)  will  be  mound-shaped  and 

approaching  normality  for  (large)  sample  sizes  such  that 

P  *  2/p (l-p7  lies  in  the  interval  (0,1)  [Ref  9:  270]. 

M 

By  applying  the  asymptotic  normality  property  of  a  binomial 
probability  distribution,  one  can  determine  the  number  of 

A 

"attacks",  M,  needed  to  have  a  confidence  interval  about  p 
of  ±.01  at  a  confidence  level  of  .99*. 


M  =  (2 . 576) 2  p(l-p) 
(•01) 2 

Yet  for  the  subproblem,  p  =  PC  thus 

M  =  (2, 576)  2  (PC) (PC) 
( • 01) 2 


(3-19) 


(3-20) 


*These  values  were  arbitrarily  chosen  and  can  be  easily  changed 
(see  Appendix  C,  Subroutine  ITER). 


Summary 

For  a  given  aim  point  strategy,  two  methods  for 
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approximating  PU  were  discussed.  First  by  using  an  event- 
composition  method,  P(E),  the  probability  of  not  accomplishing 
a  cut,  was  derived  from  the  unions  and/or  intersections  of  its 
subsets.  And  second,  by  the  use  of  a  Monte  Carlo  simulation, 
one  can  approximate  PT  by  p,  an  estimate  of  the  probability 
of  not  accomplishing  a  cut. 
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IV.  A  Search  Routine  for  Determining 
an  Aim  Point  Strategy 

Introduction 

Both  approximation  methods  from  Chapter  III,  the  event- 
composition  method  and  the  Monte  Carlo  simulation  method,  re¬ 
quire  a  given  aim  point  strategy  prior  to  estimating  the 
probability  of  not  making  a  cut,  PC.  The  second  phase  of  this 
research  is  to  develop  methods  to  determine  "good"  aim  point 
strategies.  This  chapter  reports  on  the  development  of  a 
search  algorithm  which  derives  an  aim  point  strategy  for 
achieving  the  desired  PIT.  The  algorithm  decisions  always 
choose  the  strategy  with  the  fewest  number  of  weapons,  but 
this  final  strategy  may  not  be  the  optimal  aim  point  strategy, 
i.e.  it  may  not  minimize  the  total  number  of  weapons  required. 
The  search  algorithm  is  ultimately  based  on  two  assumptions. 
Given  all  other  conditions  are  equal: 

1.  The  fewer  the  number  of  aim  points,  the  better. 

2.  Aim  point  locations  which  are  symmetric  about  the 
middle  of  a  cut  are  preferred  to  locations  which 
are  not. 

The  first  assumption  can  be  justified  in  that  there  is  less  of 
a  chance  of  incorrectly  targeting  weapons  if  there  are  fewer 
aim  points  to  work  with.  From  a  human  error  standpoint,  this 


leads  to  a  more  reliable  weapon.  The  second  assumption  is 

based  on  the  method  the  Weapon  Effectiveness  Branch  of  the 
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Analysis  Division,  Air  Force  Armament  Laboratory  currently 
uses  for  aim  point  strategies  (Ref  9).  The  assumption  is  use¬ 
ful  computationally  because  it  requires  less  computer  time. 

The  search  routine  involves  the  three  traits  of  an  aim 
point  strategy- -aim  point  locations,  number  of  aim  points,  and 
number  of  weapons  per  aim  point.  To  limit  the  search  time,  a 
feasible  range  on  the  number  of  weapons  is  determined.  Within 
this  range,  the  search  is  conducted  by  varying  one  trait  at  a 
time.  First  the  number  of  weapons  per  aim  point  is  set  to  a 
common  value  (this  weapon  allocation  is  the  aim  point  group 
designator).  Next  the  number  of  aim  points  is  fixed  (this  num¬ 
ber  identifies  an  aim  point  set  within  a  group) .  Finally,  the 
locations  are  varied  until  the  largest  PS  is  achieved.  Then 
the  number  of  aim  points  is  increased  by  one,  and  the  location 
search  is  restarted.  Once  this  search  is  completed  for  all  aim 
point  sets  within  a  group,  the  range  on  the  number  of  weapons 
and  marginal  analysis  are  used  to  determine  the  best  aim  point 
set  within  this  group. 

Determining  the  Initial  Bounds  on 
the  Required  Number  of  Weapons 

For  a  given  W  and  V  and  assuming  CEP  =  0,  symmetric 
aim  point  locations  are  determined  which  will  satisfy  the  con¬ 
straint  in  Eq  (3-1).  This  results  in  a  lower  bound,  min,  on 
the  number  of  weapons.  Then  by  relaxing  the  CEP  assumption 
such  that  CEP  equals  the  actual  weapon's  CEP  and  allocating 
one  weapon  to  each  of  these  locations,  the  probability  of 
making  the  cut  for  a  specific  aim  point  strategy,  PS,  is 
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estimated.  If  PS  >_  PC,  then  one  has  the  optimal  aim  strategy 
and  the  subproblem  is  solved.  If  PS  <  PC,  then  through  the 
use  of  marginal  analysis  and  allocating  only  one  more  weapon 
to  the  curren  total  number  of  weapons,  an  additional  weapon 
is  targeted  for  the  aim  point  which  increases  PS  the  most. 

This  weapon  allocation  is  continued  until  the  addition  of 
the  last  weapon  gives  a  PS  >  PC.  By  adding  up  the  numbers  of 
weapons  targeted  for  each  aim  point  in  this  final  aim  point 
strategy,  one  has  an  initial  upper  bound,  max,  on  the  required 
number  of  weapons.  (See  Appendix  C,  Subroutine  BOUNDS  for  the 
computer  coding  of  this  procedure.) 


Determining  the  "Best”  Aim  Point 
Locations  for  Each  Aim  Point 
Set  Within  an  Aim  Point  Group 

In  order  to  determine  the  best  locations  for  each  aim 
point  set  within  an  aim  point  group,  the  following  rules  were 
applied:  (Note  Figs  5a  and  5b  are  logic  diagrams  for  this 

entire  procedure.) 

1.  By  definition,  all  aim  points  within  any  set  will 
have  the  same  number  of  weapons  allocated  to  each  aim  point. 

2.  At  no  time  will  the  total  number  of  weapons  allo¬ 
cated  to  an  aim  point  set  exceed  the  current  upper  bound  on 
the  number  of  weapons. 

3.  For  any  initial  upper  bound  >4,  every  aim  point 
will  be  allocated  at  least  two  weapons. 

4.  For  any  initial  upper  bound  <3,  every  aim  point 
will  be  allocated  at  least  one  weapon. 


Fig  5a.  Best  Locations  for  a  Set 
Within  a  Group 
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Fig  5b.  Best  Locations  for  a  Set 
Within  a  Group 


5.  At  no  time  will  a  weapon  be  targeted  closer  than 
one  damage  radius  from  the  edge  of  the  actual  runway.  This 
rule  guarantees  that  at  least  50  percent  of  all  the  weapons 
targeted  for  the  outer  pair  of  aim  points  will  be  100  percent 
effective.  By  effective,  it  is  meant  that  the  entire  crater 
caused  by  a  weapon  will  lie  within  the  bounds  of  the  effec¬ 
tive  runway. 

6.  No  two  aim  points  will  be  located  closer  to  each 
other  than  one  damage  radius. 

7.  All  relocation  of  aim  points  will  be  in  five-foot 
increments  towards  the  center  of  the  runway.  Five  feet  was 
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chosen  because  it  is  ten  percent  of  the  smallest  V  being  in¬ 
vestigated. 


To  initialize  the  locations  of  the  aim  points  within 
a  set,  each  symmetric  pair  of  aim  points  is  located  as  far 
as  possible  from  the  center  of  cut  without  violating  Rule  5 
or  6.  (Note:  in  order  to  maintain  symmetry,  the  middle  aim 
point  for  an  odd  aim  point  set  will  always  be  located  at  the 
center  of  the  cut.)  With  these  locations,  PSq^d  is  estimated. 
Then  the  middle  pair  of  symmetric  aim  points  is  relocated 
five  feet  inboard  while  holding  all  other  aim  point  locations 
fixed.  PS  is  re-estimated,  called  If  ^NEW  <  ^OLD’ 

these  new  locations  are  the  current  "best"  locations  for  the 


aim  point  set  and  ^0LI)  is  reset  to  The  middle  pair 

is  again  moved  five  feet  inboard  and  is  estimated  again. 

This  iterative  process  will  stop  with  either  of  two  occurrences 


lm  ^NEW  >  ^OLD 
2.  Rule  6  is  violated. 


Prior  to  discussing  the  first  occurrence,  the  following 
term  must  be  defined: 


An  aim  point  subset  is  a  specific  aim  point  set  in  which 
all  symmetric  aim  point  pairs,  except  for  the  middle  pair, 
are  assigned  specific  locations.  Any  two  subsets  within 
an  aim  point  set  will  have  at  least  one  pair  of  symmetric 
aim  points,  other  than  the  middle  pair,  that  are  not  co¬ 
located.  Each  aim  point  set  will  have  at  least  one  aim 
point  subset;  both  aim  point  set  two  and  three  will  have 
the  null  set  as  their  only  aim  point  subset. 

With  the  definition  of  an  aim  point  subset,  one  can  show  that 

P5  is  a  convex  function  for  each  aim  point  subset.  Thus,  as 
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the  middle  pair  of  aim  points  are  moved  inboard,  when  FS’  starts 
to  increase,  the  "best"  locations  for  this  aim  point  subset 
have  been  passed  and  no  further  searching  is  needed  for  this 
subset.  Once  the  "best"  locations  for  a  specific  aim  point 
subset  are  found,  the  search  algorithm  moves  the  innermost 
pair  of  aim  points  (other  than  the  middle  pair) ,  which  will 
not  cause  a  violation  of  Rule  6,  five  more  feet  towards  the 
center  of  the  cut.  All  aim  points  inboard  of  the  pair  just 
moved  are  initialized  according  to  Rule  6,  and  this  new  sub¬ 
set  is  searched  for  its  "best"  location. 

Each  time  the  "best"  locations  for  a  subset  are  deter¬ 
mined,  the  FSnew  for  this  aim  point  strategy  is  compared  to 
the  current  minimum  FS0ED  for  the  entire  aim  point  set.  If 
the  aim  point  subset's  1  FSqld »  then  FSQLD  is  reset  to 

P!>NEW,  and  the  "best"  locations  for  this  aim  point  subset 
become  the  best  locations  for  the  aim  point  set.  If  F5^EW 

>  F3qLD,  t*ie  search  continues  until  one  has  the  second  occur¬ 
rence. 

The  second  occurrence  indicates  one  of  two  things. 

First,  if  Rule  6  is  violated  before  F!>NEW  >  P1>0ED,  one  has 
too  many  aim  points  and  the  middle  pair  of  aim  points  should 
be  consolidated  into  one  aim  point,  the  center  of  cut.  Second, 
if  the  locations  of  the  outermost  pair  of  aim  points  are  such 
that  when  all  inner  aim  points  are  initialized  to  be  one 
damage  radius  apart  and  the  middle  pair  of  aim  points  violates 
Rule  6,  the  search  for  the  best  locations  of  aim  points  for 
this  aim  point  set  is  complete. 
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The  search  algorithm  then  proceeds  to  the  next  aim 
point  set  within  this  aim  point  group  and  begins  to  determine 
the  best  locations  for  this  aim  point  set.  This  whole  search 
process  stops  when  Rule  2  will  be  violated  by  investigating 
the  next  aim  point  set.  (See  Appendix  C,  Subroutine  STRAT 
for  the  computer  coding  of  this  procedure.) 

Determining  the  Best  Aim 
Point  Set  Within  an  Aim 
Point  Group 

Once  the  best  locations  are  determined  for  each  aim 
point  set  (see  Fig  6) ,  the  algorithm  then  determines  which 
aim  point  set  requires  the  fewest  total  number  of  weapons  in 
order  to  have  PS  _>  PC.  This  is  accomplished  similarly  to  the 
algorithm  which  determined  the  initial  upper  bound.  Through 
marginal  analysis  an  additional  weapon  is  targeted  for  the 
aim  point  within  the  aim  point  set  which  will  cause  the 
largest  increase  in  PS.  After  each  weapon  is  allocated,  PS 
is  compared  to  PC.  Any  of  three  stopping  criteria  are  used. 
First,  if  total  number  of  weapons  allocated  is  equal  to  the 
current  upper  bound  and  PS  <  PC,  then  this  aim  point  set  is 
disregarded  and  the  next  aim  point  set  is  tested.  Second,  if 
PS  >_  PC  and  the  total  number  of  weapons  allocated  equals  the 
lower  bound  on  the  number  of  weapons  required,  the  algorithm 
stops  and  one  has  an  optimal  aim  point  strategy.  Third,  if 
PS  PC  and  the  total  number  of  weapons  allocated  in  less 
than  or  equal  to  the  current  upper  bound,  max,  this  aim  point 
strategy  and  its  PS  are  saved.  Also,  max  is  reset  to  the 
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Fig  6.  Best  Aim  Point  Set  Within 
An  Aim  Point  Group 
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number  of  weapons  required  for  this  aim  point  strategy. 

Each  aim  point  group  is  tested  and  the  one  requiring 
the  fewest  number  of  total  weapons  becomes  the  best  aim  point 
set  within  this  aim  point  group.  If  two  or  more  aim  point 
sets  tie  as  to  the  fewest  number  of  weapons,  the  aim  point 
set  with  the  highest  PS  is  then  chosen  as  the  best  aim  point 
set.  (See  Appendix  C,  Subroutine  BEST  for  the  computer  coding 
of  this  procedure.) 

The  search  algorithm  then  increases  the  number  of 
weapons  allocated  per  aim  point  by  one  and  again  starts  deter¬ 
mining  the  best  locations  for  the  aim  point  sets  within  this 
new  aim  point  group  (see  Fig  7).  This  step  is  accomplished 
because,  for  example,  the  best  locations  for  two  aim  points 
with  two  weapons  per  aim  point  may  not  be  the  best  locations 
for  two  aim  points  with  three  weapons  per  aim  point  (assuming 
that  the  current  upper  bound  is  greater  than  or  equal  to  six 
weapons).  And  in  fact,  for  the  following  input  variables: 

X  =  8000  feet 

Y  =  150  feet 

W  =  500  feet 

U  *  2000  feet 

V  *  50  feet 
CEP  =  50  feet 

P  *=  .8 

the  best  locations  for  aim  point  group  two  and  aim  point  set 
two  are  38  feet  from  the  leading  edge  of  the  runway  and  112 
feet  from  the  leading  edge  of  the  runway.  Yet  with  aim  point  group 
three  and  aim  point  set  two,  the  best  locations  are  48  feet 
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Fig  7.  The  Best  Aim  Point  Strategy 


from  the  leading  edge  of  the  runway  and  101  feet  from  the 
leading  edge  of  the  runway.  Thus,  one  cannot  assume  that  the 
"best”  locations  for  an  aim  point  set  within  a  specific  aim 
point  group  are  the  best  locations  for  that  set  within  any 
group.  (See  Appendix  C,  Subroutine  COUNT,  for  the  computer 
coding  which  determines  the  aim  point  group  and  also  stores 
the  current  best  aim  point  strategy  for  all  aim  point  groups.) 
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Methods  for  Decreasing  the 
Search  Algorithm  Calculation 
Time 

Two  methods  were  investigated  to  decrease  the  search 
algorithm's  calculation  time.  The  first  method  was  to  use 
the  quicker  of  the  two  T>(T  approximations  for  a  given  aim 
point  strategy  (either  the  event-composition  method  or  the 
Monte  Carlo  simulation) .  The  second  method  was  to  compare 


the  probability  of  damaging  the  first  takeoff  width  with  that 
of  the  middle  takeoff  width  for  a  given  aim  point  strategy. 

In  using  the  first  method  to  increase  the  efficiency 

(decrease  the  time  required)  of  the  search  algorithm,  the 

following  set  of  input  variables  were  used  to  compare  the 

computation  times  for  the  two  approximations: 

W  *  250  pounds 
CEP  =  100  feet 
X  =  4000  feet 
U  =  2000  feet 

V  =  100  feet 
P  =  .8 

10  aim  points  equally  spaced  from 
one  edge  of  the  runway  to  the  other 

two  weapons  per  aim  point 

Y  =  was  allowed  to  vary 

Since  the  driving  factor,  timewise,  for  the  event- 
composition  method  is  the  number  of  L^'s  and  requiring  that 
the  Monte  Carlo  simulation  have  2770  "attacks"  per  estimation 
(the  number  of  attacks  which  would  provide  a  confidence  inter¬ 
val  of  ±.01  at  a  confidence  level  of  .99  for  a  problem  requiring 
five  cuts) ,  it  was  assumed  that  by  varying  the  minimum  launch 
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window's  width,  one  could  determine  the  numbers  of  L.'s,  NE 
for  which  the  event-composition  method  was  always  faster  than 
the  simulation  method.  This  data  indicated  that  the  event- 
composition  method  was  always  faster  for  NE  12.  Even  though 
this  rule  (NE  _<  12;  implies  event-composition  method  will  be 
faster)  is  not  always  true,  this  fact  was  not  known  until  the 
model  was  verified. 

As  to  the  second  method  for  a  more  efficient  search 
algorithm,  given  the  "best"  locations  within  an  aim  point 
subset,  the  probability  of  damaging  the  first  minimum  launch 
window  and  the  middle  minimum  launch  window  were  assumed  to  be 
approximately  equal  to  each  other.  Thus,  by  comparing  these 
two  probabilities,  one  could  reject  some  aim  point  strategies 
from  consideration  without  having  to  use  either  method  for 
approximating  FIT.  The  author  chose  to  use  a  difference  <_  .05 
as  being  approximately  equal  and  thus  worthy  of  further  in¬ 
vestigation.  (See  Appendix  C,  Subroutine  APROX  for  the  com¬ 
puter  coding  of  this  procedure.) 

Summary 

To  use  either  of  the  two  approximating  methods  for  FU, 
one  must  first  determine  an  aim  point  strategy.  This  chapter 
discussed  a  search  routine  which  not  only  determines  aim 
point  strategies,  but  also  solves  for  the  aim  point  strategy 
which  minimizes  the  total  number  of  weapons  required  subject 
to  the  inherent  limitations  of  this  search  algorithm.  Finally, 
two  methods  for  increasing  the  efficiency  of  the  search  routine 
were  discussed. 
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V.  Verification  and  Validation 


of  the  Model 


Introduction 


Coupling  the  two  methods  for  approximating  FU  with  the 
search  algorithm,  one  has  a  model  which  will  determine  a  good 
aim  point  strategy.  Yet  this  model  is  of  questionable  value 
until  it  has  been  verified  and  validated.  Fishman  and  Kiviat 
(Ref  6)  define  verification  and  validation  as: 

Verification  -  insuring  that  the  model  behaves  the  way 
the  experimenter  intends 

Validation  -  testing  the  agreement  between  the  be¬ 
havior  of  the  model  and  that  of  the  real 
system. 

Shannon  further  expands  on  validation  by  stating,  "Validation 
is  the  process  of  bringing  to  an  acceptable  level  the  user's 
confidence  that  any  inference  about  a  system  derived  from  the 
simulation  (model)  is  correct  [Ref  12:  29]."  Furthermore,  he 
writes 

The  concept  of  validation  should  be  considered  one  of 
degree  and  not  one  of  an  either-or  notion;  it  is  not  a 
binary  decision  variable  where  the  model  is  valid  or 
invalid.  It  is  not  at  all  certain  that  it  is  ever 
theoretically  possible  to  establish  if  we  have  an 
absolutely  valid  model;  even  if  we  could,  few  managers 
would  be  willing  to  pay  the  price  [Ref  12:  208]. 

The  procedure  used  to  verify  and  validate  the  model 
was  to  first  verify  if  the  model  behaves  as  intended,  secondly 
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to  verify  the  part  of  the  model  which  estimates  PC,  and 
thirdly  to  "validate"  the  model  by  comparing  its  results  to 
the  Airfield  Assessment  Program,  AAP  (formerly  called  "MASSIVE") 
used  by  the  Weapons  Effectiveness  Branch,  Analysis  Division, 

Air  Force  Armament  Laboratory,  Eglin  AFB,  Florida  (Ref  3). 

Verification 

All  the  subroutines  within  the  model  were  verified 
through  the  use  of  print  statements  and  hand  calculations  to 
determine  whether  they  "behave"  in  the  manner  desired.  Yet, 
the  overall  results  of  the  model,  the  recommended  aim  point 
strategies,  are  extremely  dependent  upon  the  accuracy  of  the 
two  approximation  methods  for  PC.  Thus,  these  two  approxima¬ 
tion  methods  received  a  more  in-depth  verification  than  the 
rest  of  the  model. 

Verifying  the  Two  Approximating  Methods  for  PC.  After 
the  FORTRAN  coding  for  the  event-composition  method  and  the 
Monte  Carlo  simulation  method  was  shown  to  do  what  it  was 
supposed  to  do,  the  designer  checked  the  accuracy  of  both 
methods.  Since  the  event- composition  method  is  a  discrete 
approximation  of  a  continuous  process,  a  discrete  probability 
situation  was  used  in  its  initial  verification. 

Figures  8a,  b,  and  c  show  the  cases  tested  against  the 
event-composition  approximation.  In  each  case,  the  actual 
value  of  P(E),  the  probability  of  making  the  cut,  was  calculated 
manually  and  P(E),  the  probability  of  not  making  the  cut,  was 
calculated  using  the  event-composition  method.  The  results  in 


39 


Figs  8a,  b,  and  c  indicate  confidence  in  this  calculation  of 
P(E) .  One  should  note  the  following  prior  to  reviewing  the 
figures : 

1.  The  values  between  the  hash  marks  are  the  proba¬ 
bilities  that  a  weapon  targeted  for  an  aim  point  will  impact 
between  these  hash  marks. 

2.  The  number  of  hash  lines  is  indicative  of  the  num¬ 


ber  of  aim  points  (A.P.). 

3.  If  two  weapons  are  targeted  for  the  same  point, 
each  has  the  impact  probabilities  shown. 

Once  this  simplified  event- composition  approximation 
was  verified,  this  method  was  used  to  verify  the  Monte  Carlo 
simulation  approximation.  A  subproblem  was  designed  with  each 
event  composed  of  exactly  two  event  where  the  length 

over  which  an  A^^  might  occur  equals  the  standard  deviation 
of  the  weapon  being  employed.  Along  with  these  specifications, 
the  following  input  variables  were  used: 

R  =  10  feet 

Y  =  70  feet 

V  =  20  feet 

Two  aim  points,  each  located  5 
feet  either  side  of  the  center 
of  the  runway's  width 

One  weapon  targeted  for  each 
aim  point 

The  PfAj^p’s  were  calculated  using  the  formulation  discussed 
in  Chapter  III.  Both  the  event-composition  method  and  the 
Monte  Carlo  simulation  method  were  given  these  inputs,  and 
following  results  were  obtained: 
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Fig  8.  Comparison  of  P(F)  to  P(E) 


41 


P(E)  for  the  event-composition  =  .732519660803 

P(E)  for  the  simulation  =  P ( .  7  206  <_  x  £  .7434)  =  .99 
with  10,000  "attacks";  the  point  estimation  for 
the  simulation's  P(E)  =  .732 

This  seemed  to  verify  the  simulation,  but  one  more  check  was 

run.  This  next  verification  step  used  the  following  input  data 

W  =  250  pounds 
CEP  =  50  feet 
X  =  8000  feet 

Y  =  150  feet 
U  =  2000  feet 

V  =  100  feet 
P  =  .8 

The  entire  model  was  run  31  times  using  a  uniform  random  num¬ 
ber  generator  to  provide  a  different  "seed"  for  the  normal 
random  number  generator  on  each  run.  This  verification  step 
was  taken  in  order  to  demonstrate  consistency  in  the  simulation 
method.  The  following  results  were  obtained: 

the  mean  of  the  31  observations,  X  =  .960429 
the  unbiased  sample  variance,  S  =  .00505 

Assuming  the  Central  Limit  Theorem  applied  to  these  results 
(Ref  11:  252),  one  would  expect  about  21  observations  within 
one  standard  deviation  of  X  and  about  29  observations  within 
two  standard  deviations  of  X.  The  actual  results  had  25  obser¬ 
vations  within  the  first  interval  and  all  31  observations 
within  the  second  interval.  Due  to  the  smallness  of  the  number 
of  observations  taken,  these  discrepancies  were  not  considered 
severe,  and  the  simulation  method  was  assumed  to  be  consistent 
and  not  "seed"  dependent. 
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With  the  verification  of  the  simulation  method  completed, 
a  number  of  test  cases  were  run  to  compare  the  simulation  method 
and  the  event-composition  method.  (Note:  when  the  simulation 
method  was  verified  using  the  event- composition  method,  the 
subproblem  was  constructed  using  a  very  simplified  event- 
composition.  That  simplification  was  not  based  on  the  same 
size  Z's  described  in  Chapter  III.  However,  this  comparison 
will  actually  compare  the  simulation  and  the  event-composition 
methods  as  developed  in  Chapter  III.)  All  input  variables  were 
held  constant  except  for  W,  CEP  and  Y.  Furthermore,  in  each 
case  there  were  four  aim  points  located  at  .2Y,  .4Y,  . 6Y  and 
. 8Y ;  and  one  weapon  was  targeted  for  each  aim  point  (see  Table 

I). 

TABLE  I 


Comparison  of  PIT  Approximation  Methods 


Inputs 

PtT  Approximations 

W 

CEP 

Y 

Event-Composition  Simulation 

f . 9 5  confidence  interval 

1000 

20 

300 

. 000484 

70 

to  .061980 

1000 

100 

150 

.051687 

.061378 

to  .094621 

1000 

100 

300 

.388422 

. 411219 

to  .472781 

250 

20 

300 

.001631 

.0031465 

to  .014853 

X  * 

4000, 

U  =  2000,  V  =  100 

The  results  of  these  comparisons  indicated: 

1.  Though  the  event-composition  method  was  designed 
to  be  "optimistic",  i.e.  give  a  larger  PU  than  the  actual 
attack,  this  approximation  seems  to  be  "pessimistic." 

2.  Within  the  range  of  PC  under  investigation,  .946 
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to  .997  (P  =  .8;  CEP's  such  that  one  has  four  or  five  cuts), 
the  event- composition  method  would  probably  not  be  enough  in 
error  to  suggest  using  fewer  weapons  than  would  actually  be 
required.  (This  assumption  was  shown  to  be  true  for  the  range 
of  input  variables  investigated.  See  Table  II) 

The  conclusion  derived  from  these  results  was  to  use 
the  event-composition  method  in  the  search  algorithm  when  it 
provided  a  time  savings,  but  to  always  use  the  Monte  Carlo 
simulation  method  to  approximate  the  resulting  PC  from  the 
best  aim  point  strategy.  This  procedure  would  provide  a 
built-in  model  verification  for  any  aim  point  strategy  deter¬ 
mined  "best"  by  the  event-composition  method. 

Verification  of  the  Time  Saving  Methods.  As  was  indi¬ 
cated  in  Chapter  IV,  using  12  I^’s  as  the  changeover  point 
between  the  efficiency,  timewise,  of  the  event-composition 
method  and  Monte  Carlo  simulation  turned  out  to  be  in  error. 
Using  the  following  set  of  input  data: 

W  =  1000  pounds 
CEP  =  20  feet 
X  =  8000  feet 

Y  =  150  feet 
U  =  2000  feet 

V  =  60  feet 
P  =  .8 

3  Symmetric  Aim  Points 
2  Weapons  per  Aim  Point 
3470  "Attacks"  for  the  simulation 

the  simulation  method  was  ten  times  faster  than  the  event- 

composition  method.  Due  to  lack  of  time,  the  relationship 
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between  the  number  of  weapons  (the  major  determinant  of  the 
time  for  the  simulation)  and  the  number  of  L^'s  (the  driving 
factor,  timewise,  for  the  event- composition  method)  was  not 
investigated  further. 

The  second  efficiency  method  used  in  the  search  routine 

setting  a  tolerance  on  the  difference  between  the  probability 

of  destroying  the  first  V  and  the  probability  of  destroying 

the  middle  V,  was  verified  with  the  following  data: 

W  =  500  pounds 
CEP  =  50  feet 
X  =  8000  feet 

Y  =  150  feet 
U  =  2000  feet 

Y  =  50  feet 
P  =  .8 

The  above  data  was  run  with  two  tolerances,  the  actual  toler¬ 
ance  used  in  the  model,  .05,  and  the  tolerance  of  1.0.  Both 
of  these  computer  runs  determine  the  same  best  aim  point 
strategy,  yet  the  run  with  a  tolerance  of  .05  was  1.4  seconds 
(approximately  one  percent)  quicker. 

3 

Factorial  Design.  Due  to  the  lack  of  time,  only  a  2 
factorial  design  was  constructed  to  investigate  the  model. 

The  following  input  variables  were  held  constant  for  each  run 
of  the  model: 

X  =  8000  feet 

Y  =  150  feet 
U  =  2000  feet 
P  =  .8 

The  remaining  input  variables  were  set  at  the  following  levels 


45 


Low 


High 

W  1000  pounds  250  pounds 

CEP  100  feet  20  feet 

V  100  feet  50  feet 

Table  II  provides  the  levels  of  the  varied  inputs  and  the 
model's  results  for  each  run.  Using  the  number  of  weapons 
required  to  make  the  cut  as  the  dependent  variable  and  assign¬ 
ing  qualitative  value  (zero  for  low  and  one  for  high)  to  each 
of  the  three  independent  variables,  an  analysis  of  variance 
(ANOVA) ,  was  performed  "...  to  locate  important  independent 
variables  in  a  study  and  determine  how  they  interact  and 
affect  the  response  (the  dependent  variable)  [Ref  11:458]." 

In  order  to  accomplish  the  ANOVA,  the  three-way  interaction 
of  the  independent  variables  had  to  be  suppressed.  This  sup¬ 
pression  was  caused  by  too  little  data  to  investigate  the 
main  effects  and  all  their  interaction.  The  result  of  this 
ANOVA  suggested  the  only  significant  effects  on  the  number 
of  weapons  required  were  caused  by  W,  CEP,  and  V  and  not  their 
interactions.  Thus,  a  second  ANOVA  was  accomplished  suppres¬ 
sing  all  interactions.  The  result  of  this  ANOVA  are  given 
in  Table  III.  Table  III  indicates,  for  a  Type  I  error  (Ref  11 
328)  of  less  than  .05,  the  only  independent  variables  (of 
those  investigated)  which  definitely  affect  the  number  of 
weapons  required  are  CEP  and  V. 

With  the  insight  provided  by  the  ANOVA  in  Table  III, 
a  second  type  of  analysis  was  attempted.  This  analysis  was 
a  linear  regression  using  the  same  dependent  and  independent 


Results*  from  the  2  Factorial  Design 


the  estimated  probability  of  destroying  the 


TABLE  III 

ANOVA  of  the  Main  Effects  of  W,  CEP, 
and  V  on  the  Number  of  Weapons 
Required  to  Make  a  Cut 


Source  of 
Variation 

Sum  of 
Squares 

Degrees 

of 

Freedom 

Mean 

Square 

Error 

Significance 

F  of  F 

W 

3.125 

1 

3.125 

2.273 

.206 

CEP 

45.125 

1 

45.125 

32.818 

.005 

V 

28.125 

1 

28.125 

20.455 

.011 

Residual 

5.500 

4 

1.  375 

Total 

81.875 

7 

11.696 

with  X  =  8000,  Y  = 

150,  U  = 

2000,  and 

P  =  .8 

variable.  Yet,  the  independent  variables  were  now  treated  as 
quantitative  factors  rather  than  qualitative,  as  in  the  ANOVA. 
The  following  linear  equation  was  obtained: 

the  number  of  weapons  required  = 

. 059375CEP  -  .075V  +  6.6875  (5-1) 

The  coefficients  of  CEP  and  V  and  the  constant  term  are  all 
significant  for  a  Type  I  error  <_  .01  (that  is,  one  is  98  per¬ 
cent  confident  that  none  of  these  values  are  actually  zero) , 

2 

and  for  the  adjusted  R  value  89.5  percent  of  the  variance  in 
the  number  of  weapons  required  is  explained  by  Eq  (5-1). 

"Validation"  of  the  Model 

As  was  pointed  out  in  the  introduction  to  this  chapter, 
it  may  very  well  be  impossible  to  establish  a  model  as  absol¬ 
utely  valid  (Ref  12:  208).  Yet,  Shannon  states,  ”...  it  is 
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the  operational  utility  of  the  model  and  not  the  truth  of  its 
structure  that  usually  concerns  us  [Ref  12:  29]."  It  was  felt 
that  the  best  measure  of  the  "operational  utility"  of  this 
model  would  be  to  compare  its  results  with  a  model  such  as 
AAP  which  is  currently  being  used  by  the  Air  Force  for  similar 
types  of  investigation  (Ref  3) . 

Before  discussing  the  comparison,  one  should  understand 
that  AAP  was  designed  to  determine  the  probability  of  destroy¬ 
ing  a  runway  given  the  aim  point  strategy.  During  conversa¬ 
tions  with  Mr.  Bass  (Ref  3) ,  it  became  evident  that  to  obtain 
similar  results,  he  would  need  to  run  a  number  of  different 
aim  point  strategies.  Plus,  the  aim  point  strategies  he  em¬ 
ployed  would  have  to  be  based  on  his  intuition  gained  from 
dealing  with  similar  problems  in  this  field.  (Note:  for  an 
in-depth  discussion  of  AAP,  see  Ref  2.) 

As  to  the  comparison  of  the  results  of  the  two  models, 

the  following  input  data  were  used: 

W  =  500  pounds 
CEP  =  50  feet 
X  =  8000  feet 

Y  =  150  feet 
U  =  2000  feet 

V  =  50  feet 
P  =  .8 


Bass  ’ 

"best" 

aim  point  strategy  was: 

1. 

Total 

#  of  weapons  was  24 

2. 

Three 

weapons  were  on  each  of 

two  aim 

points 

3. 

These 

aim  points  were  located 

50  feet 

from  one 
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edge  of  the  runway  and  50  feet  from  the  other  edge. 

4.  Four  cuts  were  required  to  be  made  in  the  runway. 

5.  The  cuts  were  located  at  1600-foot  intervals  along 
the  runway  (1600,  3200,  4800,  6400  feet  from  the 
end  of  the  runway) . 

6.  The  estimated  probability  of  destroying  the  runway 
was  .88. 

This  model  obtained  the  following  results: 

1.  Total  #  of  weapons  was  24 

2.  Three  weapons  on  each  of  two  aim  points 

3.  These  aim  points  were  located  48  feet  from  one  edge 
of  the  runway  and  49  feet  from  the  other  edge  (lack 
of  symmetry  is  due  to  rounding  the  aim  point  loca¬ 
tions  off  to  whole  numbers  of  feet) 

4.  Four  cuts  were  required  in  the  runway 

5.  The  cuts  were  located  at  a  distance  of  1865  feet, 
3605  feet,  5345  feet,  and  7084  feet  from  the  end 
of  the  runway 

6.  The  point  estimate  for  the  probability  of  destroy¬ 
ing  the  runway  was  .908 

This  aim  point  strategy  was  then  employed  using  AAP,  and  the 
point  estimate  for  the  probability  of  destroying  the  runway 
increased  to  .89  with  a  standard  deviation  ot  .314  (only  100 
"attacks"  were  conducted) . 

Despite  the  fact  that  only  one  comparison  was  made 
between  the  model  developed  in  this  thesis  and  AAP,  it  seems 
clear  that  there  is  some  "operational  utility"  in  the  model. 
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This  comparison  showed  this  model  determined  a  "better"  aim 
point  strategy  than  the  one  derived  by  Mr.  Bass'  intuition 
and  AAP.  Yet,  one  comparison  cannot  be  used  as  validation 
for  this  model.  What  can  be  ascertained  is  that,  for  the 
case  investigated,  the  results  of  the  model  are  not  invalid. 
Further  cases  should  be  tested  before  one  could  suggest  that 
the  model  is  valid,  but  due  to  time  constraints  on  the  AFIT 
computer  and  Mr.  Bass'  schedule,  the  investigation  was  not 
continued. 

Summary 

This  chapter  discussed  the  verification  and  validation 
of  the  model.  The  steps  taken  in  verifying  the  model  were  1) 
making  sure  that  each  of  the  subroutines  accomplished  what  the 
designer  intended  it  to;  2)  an  in-depth  verification  of  both 
methods  of  approximating  PTT,  given  an  aim  point  strategy; 
and  3)  using  an  ANOVA  and  linear  regression  to  determine  the 
effects  of  W,  CEP  and  V  on  the  number  of  weapons  required  to 
achieve  PC.  The  "validation"  consisted  of  comparing  the 
results  of  this  model  to  those  of  AAP  for  a  specific  set  of 
input  variables. 
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VI.  Conclusions  and  Recommendations 


This  study  reports  the  development  of  a  model  which 
determines  an  aim  point  strategy  for  reducing  the  total 
number  of  weapons  employed  to  achieve  a  desired  level  of 
runway  destruction.  Initially,  the  overall  problem  is  re¬ 
duced  to  a  number  of  identical,  stochastically  independent 
subproblems.  Using  one  of  these  subproblems,  two  methods 
for  estimating  its  destruction  criteria  are  developed. 

Then,  these  estimating  methods  are  coupled  with  a  search  algor¬ 
ithm  to  determine  the  "best"  aim  point  strategy.  Finally,  the 
entire  model  is  verified  and  "validated."  From  the  effort 
expended  on  this  study,  a  number  of  conclusions  and  recom¬ 
mendations  can  be  made. 

Conclusions 

1.  The  model,  when  compared  to  the  current  Air  Force 
work  in  this  area,  seems  to  provide  aim  point  strategies  which 
are  better  than  those  currently  used. 

2.  The  results  of  this  study  indicate  an  area  of  pos¬ 
sible  savings  for  the  Air  Force.  This  savings  is  brought 
about  by  reducing  the  number  of  weapons  employed  while  main¬ 
taining  the  desired  level  of  destruction. 

3.  This  study,  in  its  entirety,  provides  a  basis  for 
future  investigation  on  targeting  conventional  weapons  against 
ron-reinforced  concrete  runways. 
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Recommendations 


1.  Since  symmetric  aim  points  are  currently  accept¬ 
able  in  determining  an  aim  point  strategy  (Ref  3) ,  the 
approximation  methods  could  be  changed  to  incorporate  this 
type  of  strategy  and,  thus,  further  increase  their  efficiency 
timewise . 

2.  Further  comparisons  between  AAP  and  this  model 
need  to  be  performed  to  increase  the  user's  confidence  level 
in  the  model's  results. 

3.  One  could  adapt  the  model  to  investigate  CEPs 
which  cannot  logically  be  described  by  circular  normal  prob¬ 
ability  distributions. 

4.  One  could  adapt  the  model  to  use  different  aim 
point  strategies  for  some  cuts  and  possibly  further  reduce 
the  total  number  of  weapons  required  to  achieve  the  runway 
destruction  criterion. 

5.  Due  to  the  availability  of  data,  this  model 
determines  aim  point  strategies  for  destroying  non-reinforced 
concrete  runways.  A  slight  modification  in  Subroutine  CONVRT 
would  allow  investigations  of  reinforced  concrete  runways 
which  are  more  typical  of  current  operational  runways. 
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Derivation  of  the  Effective  Explosive 


Weight  Conversion  Formula 

For  a  known  effective  explosive  weight,  W,  the  optimal 
DOB  can  be  determined  from  the  following  inequality  (Ref  10: 
1-133)  : 

1  <  D0B/W1/3  <  1.5  (A- 1) 


for  W  >  50  pounds  of  TNT.  Assuming  the  midpoint  of  this 

interval  (1.25)  to  be  a  good  approximation  of  this  ratio,  one 

1/3 

can  then  determine  the  scaled  radius  (ft/lb  )  by  referring 
to  Figs  5.9  and  5.10  in  Waterways  Experiment  Station  Techni¬ 
cal  Report  Number  N-76-2  (Ref  13 :100- 01) .  For  a  W  greater 
than  50  pounds  of  TNT  detonating  in  dry-to-moist  clay  beneath 

a  non-reinforced  concrete  runway  (Ref  13:100),  the  scaled 

1/3 

radius  is  approximately  3.65  ft/lb  .  For  a  W  greater  than 
50  pounds  of  TNT  detonating  in  dry-to-noist  sand  beneath  a 
non-reinforced  concrete  runway  (Ref  13:101),  the  scaled  radius 
is  approximately  3.54  ft/lb^3. 

These  scaling  radii  are  derived  from  a  scaling  factor 
determined  by  Chabai  (Ref  4) : 


p.  1/3  W  1/3 

(^)  (|T> 

p2  2 


(A-  2) 


where  R  is  the  damage  radius,  and  p  is  the  density  of  the 
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material  in  which  the  explosion  occurs.  The  subscript  1  refers 
to  a  set  of  known  data.  For  this  use  P]/p2  =  thus 

R,  W1  1/3 

n~  =  (ra— )  (A- 3) 

R2  2 

and 

R.  R- 

- 173  =  - 173  (A-4) 

(w^173  (w2)1/3 

The  expression  on  the  righthand  side  of  Eq  (A-4)  is  the  scaled 
radius.  Since  the  density  radio  (p^/p2)  is  assumed  to  be  one, 
the  scaled  radius  will  be  assumed  to  be  the  lesser  of  the  two 
from  the  previous  paragraph  (3.54  ft/lb1^3).  Thus 

R  =  (3.54)(W)1/3  (A- 5) 

For  example,  with  an  effective  explosive  weight  of 
1000  pounds  TNT,  the  damage  radius  would  be  35.4  feet. 
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APPENDIX  B 

Derivation  of  the  Circular  Error 
Probable  Conversion  Formula 
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Derivation  of  the  Circular  Error 


Probable  Conversion  Formula 

For  a  bivariate  normal  probability  distribution, 
where  the  x  and  y  variables  are  independent,  the  joint  den¬ 
sity  function  is 

f(x,y)  =  {(_l—)exp-icpi)2}{(— )exp-i(pi)2}  (B-l) 

oy/2?  L  °x  a  /?¥  L  °y 

x  y 

The  first  assumption  in  Chapter  I  is  that  the  CEP  of  the 
weapon  can  logically  be  described  by  the  circular  normal 
distribution.  This  implies  that  the  standard  deviations  of 
the  joint  distribution  are  "...  independent  of  direction 
or  choice  of  coordinate  axis.  .  ."  (Ref  4:14),  or  more  spe¬ 
cifically,  o=a.  Thus,  applying  this  assumption  to  Eq 
x  y 

(B-l),  one  obtains 

f(x,y)  =  — lI{exp-|(^)2}{exp-|(^)2}  (B-2) 

2-no 

A  further  simplification  is  achieved  by  realizing  that  y  is 
actually  the  aim  point  of  the  weapon  and,  thus,  y  =  0. 

2  2 

f(x,y)  =  — K  exp-|(£— £-)  (B-3) 

2tt  a  o 

Letting  r2  =  x2  +  y2,  Eq  (B-3)  becomes  (Ref  4:14): 

f(r)  =  .1  exp  -  |(I)2  (B-4) 

2tto 
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Yet  CEP  is  defined  such  that 


Thus , 


.5  =  /0CEP  f(r)dr 


(B-5) 


c  _  fCEP  1  .  l,r,2. 

•  5  •  /a  2  *2$  o 


'o 


2ira 


_  n  1,CEP.2 

.5  =  1-  exp  -  ) 

-  _  l,CEPs2 

2  -  exp  2(  a  ) 

.  -  _  1  ,-CEP.  2 

in  2  -  jC— ) 

2  In  2  =  C~) 2 


a  = 


CEP 


(2  In  2) 


172 


(B-6) 

(B-7) 
(B-8) 
(B-9) 
(B- 10) 
(B- 11) 
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FORTRAN 


APPENDIX  C 
Coding  for  the  Model 
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FORTRAN  Coding  for  the  Model 


PROGRAM  MAIN( INPUT, OUTPUT, TAPE5=  INPUT, TAPE6=0UTPUT) 


C 


THIS  PROGRAM  IS  OESISNEO  AROUNO  A  MODULAR  APPROACH 
HOST  OF  THE  CALCULATIONS  ARE  DONE  IN  THE  SUBROUTINES 
THE  FOLLOWING  VARIABLES  ARE* 

1.  W=THE  EFFECTIVE  EXPLOSIVE  WEIGHT  OF  THE  WEAPON 

2.  CFP=  THE  CIRCULAR  ERROR  PROBABLE  OF  THE  WEAPON 

3.  RUNLEN=THE  ACTUAL  RUNWAY  LENGTH 

4.  RUNWIOsTHE  ACTUAL  RUNWAY  WIDTH 

5.  TOLENsTHE  MINIMUM  TAKEOFF  LENGTH 

6.  TOWIO=THE  MINIMUM  TAKEOFF  WIDTH 

7.  PRORUN=THE  RUNWAY  DESTRUCTION  CRITERIA'S  CONFIDENCE  LEVEL 


C 


DIMENSION  MINAP(20)  ,MINBMB(20) 

DIMENSION  A IMLOC ( 20) , NUMB OM( 20 ) «  SECBON (20,20) » PRO SEC (20,20) 
DIMENSION  PR OF AC C20, 201, DA ML OC (20) 

DIMENSION  INDX(20),T(20),SUM(20),  SECFRAC  20,20) 

DIMENSION  LOCCUT(20)  ,  APMIN<20)  ,MINBOM(  20),  BEST  API  20,20)  ,TMRGIN<  20) 
DIMENSIO)  VALUE  (20)  ,AI  KJP  (  20, 20) 

DIMENSION  B  VALUE ( 20) 

-  DIMENSION  HITC50) ,R1( l)  * .  . . .  " 

DIMENSION  INO(20) 

DIMENSION  EQUPR0(2,  20,25)  ,S1(25)  ,  PROD(  2)  ,  PROEQUC  2, 20, 25) 

COMMON  PROSPG (12,20,20) 

- COMMON /JOHN /SECBON, PROSEC, PROFAC ,DAHLOC» SECFRA  - - 

COMMON /CUT/LOCCUT 

-  COMMON/M BR3IN/TMRG IN  - - -  -  - 

COMMON/A IMGRP/AIMGP 

"*  COMMON/I  NOICE/INOX  . . . 

COMMON/FIXE  D/RUNWID,  OAMRAD,  STADEV,  TOWID,  PROCUT 

. COMMON/OPT/BESTAP, IND,BVALUE 

COMMON/S  IHS/H  IT,  Ri 
COMM3N/ITER AT/NITERA 

COMMON/GRP/NSPG,  NSPW,NGPW,  STEP,  FRA  I,  ERUNMI,  ERWID 
COMMON/E  QUAL/MINAP,MINNBR,  MI NBMB  , NUMBER 
COMMON/ A  PPROX/EQUPRO, S i,PROO,PROEQU 


500  FORMAT  <  1H1I 

RANSET  =585285124. 

100  RE  AO*  ,  H  ,  CEP  ,  RUNLEN  ,  RUN  HD ,  TOLEN, TOWID,  PRORUN 
IFCE0FC5)  .NE.O)  GO  TO  20 


C  INPUTS  T3  THE  SYSTEM 

MRITE  <  6»  500) 

PRINT**"  ~ 

PRINT*,"  THE  EFFECTIVE  EXPLOSIVE  HEIGHT  OF  THE  WEAPON  IS  ",H,"  POU 
CNOS  OF  TNT." 

PRINT*,"  " 

PRINT*,"  THE  CEP  OF  TH E  WEAPON  IS  ",CEP,"  FEET." 

PRINT * ,"  " 

PRINT*,"  THE  ACTUAL  RUNWAY  DIMENSIONS  ARE  ", RUNLEN,"  FEET  BY  ",R 
CUNWIO, "  FEET." 

PRINT*,"  " 

PRINT*,"  THE  MINIMUM  LAUNCH  WINDOW  DIMENSIONS  ARE  ",TOLEN,"  FEET 
CBY  ",  TOWID,"  FEET." 

PRINT*,"  " 

PRINT*,"  THE  RUNWAY  DESTRUCTION  CRITERIA  IS  ", PRORUN,  "  ." 

PRINT*,"  " 

PRINT*,"  " 

C  THE  CALLS  TO  THE  SUBROUTINES  WHERE  MOST  OF  THE  WORK  IS  DONE 

CALL  CON VRT  <  W,CEP, RUNLEN, TOLEN,  PRORUN,  NCUT) 

CALL  GROUP 

CALL  ITER  _ 

. CALL  BOUNOSC  MIN, MAX, NU MAIM) 

IFCMAX.EQ.MIN)  NUM3ER=MJMAIM 

IFCMAX .EQ.MIN)  GO  TO  15 

ICOUNT=l 

MINNBR=MIN 

MINNUM*M AX 

IFUMAX.GE.3)  .AND. (MAX .LE .4) >  ICOUNT*0 
CALL  COUNT!  ICOUNT  ,  MINNUM) 
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C  THE  PRINT  STATEMENTS  FOR  THE  OUTPUT 

15  CONTINUE 

PRINT*,**  - 
MINN3R=0 

DO  21  1*1, NUMBER 

PRINT*,-  -,HIN8MBtI>  ,"  WEAPONS  SHOULO  BE  TARGE TEO  FOR  ",  NINAPI  I)  , 
C"  FEET  PROM  THE  EDGE  OF  THE  RUNWAY." 

AIMLOC<II*MINAP<I) 

MINNBRSMINNBR*MINBMB( I) 

21  CONTINUE 
PRINT*,"  - 

DO  22  1=1 ,NCUT 

PRINT*,"  THE  LOCATION  OF  CUT  NUMBER  ",I,"  IS  ",LOCCUTf I)  , "  FEET  FR 
COM  THE  LEAOING  EDGE  OF  THE  RUNWAY. " 

22  CONTINUE 
NTOTAL*NCUT*MINNBR 
PRINT*,"  - 

CALL  S IM( NUMBER, A IML0C«MI NBMB,  TOTAL! 

Y*l.-TOTAL 
Y1*Y- • 01 
Y2=Y* • 01 

IF (Y2. GT .1. !  Y2-1 • 

PRINT*,"  WITH  THIS  AIM  POINT  STRATEGY,  ONE  IS  99  PERCENT  CONFIDENT 
C  THAT  THE  ACTUAL  CUT  DESTRUCTION  CRITERIA  IS  BETWEEN  -,Y1,"  AND  -, 
CY2 

PRINT*,"  " 

PRINT*,"  THE  TOTAL  NUMBER  OF  HEAPONS  REQUIRED  IS  "»NTOTAL 
PRINT*,"  " 

’  Xi=NCUT 

CALL  RUNPRO(Y,NCUT,BNOLOW,BNOUP! 

*  IFIBNOUP  .GT.l.l  BNOUP=l. 

"  PRINT*,"  WITH  THIS  AIM  POINT  STRATEGY,  ONE  IS  99  PERCENT  CONFIDENT 
C  THAT  THE  ACTUAL  RUNWAY  DESTRUCTION  CRITERIA  IS  BET WEEN", BNOLOH, 

C"  ANO  ", BNOUP 
1000  CONTINUE 


GO  TO  100 

20  STOP  "END  OF  PROGRAM- 
END 
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Subroutine  CONVRT 


SUBROUTINE  CONVRT C W, CEP, RUNLEN*TOLEN, PRORUN, NCUTI 


THIS  SUBROUTINE  <CO*VRT)  DOES  THE  INITIAL  CONVERSION  ON  THE  INPUT 
VARIABLES.  IT  RETURNS  TO  THE  MAIN  PROGRAM  ALL  THE  INFORMATION 
REQUIRED  TO  CONTINUE  PLUS  THE  LOCATIONS  OF  THE  CUTS  ARE  ROUNDED 
OFF  TO  THE  NEAREST  INTEGER  NUMBER  OF  FEET  AND  STOREO  IN  LOCCUT 
THE  FOLLHING  VARIABLES  ARE  I 

1.  OAMRAO=THE  OAMAGE  RADIUS 

2.  STADEV=THE  STANDARD  DEVIATION 

3.  TOM=MINIMUM  TAKEOFF  LENGTH  FOR  ACROSS  THE  RUNWAY  TAKEOFFS 

4.  NCUT=THE  NUMBER  OF  CUTS  REQUIRED  IN  THE  RUNWAY 

5.  PROCUT=THE  REQUIREO  PROBABILITY  FOR  ANY  ONE  CUT 

6.  LOCCUT-AN  ARRAY  OF  THE  LOCATIONS  OF  THE  CUTS 


C 


DIMENSION  LOCCUT(  20) 

COMMON/CUT/LOCCUT 

COMMON/FIXE D/RUNWIC,OAMRAO,STAOEV» TO WTO, PROCUT 
CONNON/GRP/NSPG»NSPW, NGPW,  STEP*  FRAI, ERUNWI.ERWIO 


C 


DAMRA0  =  3.54M(W>*M1./3.I  ) 

STADEV  =  <CEP)/  I  <2.*AL0G(2.))*M  1./2.)  ) 
ERUNWI  =  RUNWIOf2.*OAMRA  f 

- TOM=< (TOLEN) **2.- C  KUNWIO) »*2.) **f 1./2.I 

BCI*T0M-6.*STA0EV 

- NCUTs  CRUNLEN- (TOM-3.  *STADEV)  1/BC  1*1 

X»NCUT 

PROCUT*  (PRORUNJ  **  *1. /X  ) 

LOCCUT  Cl  l=TOM-3.*STAOEV 
Y*LOCCUT  ( 11 
DO  1  1=2,  NCUT 
Z*I-l 

LOCCUTtII=Y*Z*BCI 
l  CONTINUE 

RETURN 
END 
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Subroutine  GROUP 


SUBROUTINE  GROUP 


C 


THIS  SUBROUTINE  (GROUP)  IS  OESIGNEO  TO  PROVIOE 
THE  NECESSARY  INFORMATION  FOR  THE  EVENT  COMPO¬ 
SITION  METHOD  AND  ALSO  TO  COMPUTE  THE  NUMBER 
OF  GROUPS  PER  WIDTH  (NGPW). 

THE  FOLLOWING  VARIABLES  ARE  I 

1.  STEP*THE  SIZE  OF  TH f  SECTION 

2.  ERUNWI=THE  EFFECTIVE  RUNWAY  WIDTH 

3.  ERWID=TH£  EFFECTIVE  RUNWAY  WIDTH  MINUS  THE 
LAST  CAMAGE  RADIUS 

4.  NSPH=NUMBER  OF  SECTION  PER  WIDTH 

5.  NSPG=NUMBER  OF  SECTIONS  PER  GROUP 

6.  FRAC=THE  FRACTION  OF  TOWIO  WHEN  COMPARED  TO 
THE  DAMAGE  RADIUS 

7.  FRAsTHE  FRACTION  OF  THE  GROUP  WHEN  COMPARED 
TO  STEP 

9.  NGPW=  NUMBER  OF  GROUPS  PER  WIDTH 


COMMON/FIXED/RUNW 10, DAMRAO,STAOEV» TOWIO, PROCUT 
C0HM0N7GRP/ NSPG, NSPW, N  CPW , STEP ,FRAI, ERUNWI, ERWID 


W5,#1  -  - 
STEP*TOWIO*A85 
FRAC* TOW ID/ OAMRAO 
NFRAC’TOWlO/OAMRAO 
'  Y*NFRAC 

FRACI=FRAC-Y 
OAHFRA*FRACI*OAMRAO 
PRINT*,-  FRACI=  - , FRAC  I 
PRINT*,"  OAHFRAs  ", OAM FRA 

C  EFFECTIVE  RUNWAY  WIDTH  (ERUNWI) 

ERUNWI>RL'NWI0*2.*0AMRAC 
ERWID*  RUNWI D ♦ OAMRAO 
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C  NUMBER  OF  SECTIONS  PER  GRO  tP  CNSPG) 

NSPG*  (  OA hRA  D^TOMID^OAH  FRA  I /STEP* 1 
C  NUMBER  OF  SECTIONS  PFR  WIOIH  INSPN) 
NSPM*ERN  IO/STEP+1 

C  NUMBER  OF  GROUPS  PER  MIOTH  (NGPM) 
NGPH=N  SPW-NSPG+1 

C  WORKING  WITH  THE  “FRACTION" 

FRA*<OAMRAO*TOWID*OAMFRAI  /STEP 
NFRAsCDAMRAO»TOMIO*OAMFRAI/STEP 
V*NFR A 
FRAI*FRA-Y 

return 

ENO 
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Subroutine  ITER 


f 

lr 


SUBROUTINE  ITER 


THIS  SUBROUT  INE  <  ITER)  CALCULATES  THE  NUMBER  OF 
ITERATIONS  (NITERA)  WHICH  ARE  REQUIRED  BY  THE  SIMULATION 
FOR  ONE  TO  BE  <m  CONFIDENT  THAT  THE  ACTUAL  PROBABILITY  OF 
MAKING  THE  CUT  IS  WITHIN  .Cl  OF  THE  VALUE  COMPUTED 


C 


COMMON/F IXE O/RUNW I D» DAMRA  D»STADE  V»TOWIO»  PROCUT 
"  COMMON/ITERAT/NITERA 


X*  C 12.  576**2.)  *<  PROCUT  >*<  1.- -PROCUT)  !/«(•  01)  *  *2. » 
NITERA*X  *1. 

RETURN 

END 
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Subroutine  BOUNDS 


SUBROUTINE  BOUNDS < MIN, MAX  , NUMA IM> 


C 


THIS  SUBROUTINE  (BOUNOSI  COMPUTES  THE  MINIMUM  LONER 
BOUNOS  ON  THE  WEAPON  AND  AN  INITIAL  UPPER  BOUND,  THESE 
BOUNDS  ARE  ON  THE  NUMBER  OF  WEAPONS  REQUIRED  TO  MAKE 
A  CUT  WITH  THE  GIVEN  DESTRUCTION  CRITERIA  UN  THIS 
CASE  PROCUTI 

THE  FOLLOWING  VARIABLES  ARE  I 

1.  MINsTHE  MINIMUM  NUMCER  OF  WEAPONS  REQUIRED  (THIS 
NUMBER  COMES  FROM  USING  AN  AIM  POINT  STRATEGY 
WHICH  WILL  GUARANTEE  THAT  NO  SECTION  OF  LENGTH  TOWID 
WILL  REMAIN  UNDAMAGED.  THEN  MIN  IS  THE  NUMBER 

AIM  POINTS  FOR  THE  STRATEGY  CONE  IS  ASSUMING  THAT 
THE  CEP  IS  0) 

2.  MAX=THE  NUMBER  OF  WEAPONS  REQUIRED  IF  THE  MIN 
AIM  POINT  STRATEGY  IS  USED  WITH  THE  ACTUAL  CEP 
OF  THE  WEAPON  UNOER  INVESTIGATION 

3.  TMARGIN  IS  AN  ARRAY  FOR  USE  IN  COMPARING  THE 
MARGINAL  IMPROVEMENT  IN  THE  PROBABILIT YY  OF 
MAKING  THE  CUT  WITH  ONE  MORE  WEAPON  ADDED  TO 
EACH  AIM  POINT.  THIS  IS  DONE  SEPARATELY  SO  THAT 
ON  NO  ONE  ITERATION  WILL  THE  TOTAL  NUMBER 

OF  WEAPONS  INCREASE  BY  MORE  THEN  ONE 
%.  PROMIN=THE  PROBABILITY  OF  MAKING  THE  CUT 
WITH  THE  MIN  AIM  POINT  STRATEGY 

5.  AIMLOC=AN  ARRAY  OF  THE  AIM  POINT  LOCATIONS 

6.  EPS ILN=A  CONVERGENCE  CRITERIA  FOR  SUBROUTINE  UNION 
WHICH  IS  SET  TO  ZERC  FOR  THE  USE  IN  THIS  THESIS 

7.  NJMAI  M=THE  NUMBER  OF  AIM  POINTS 


DIMENSION  A I MLOC (  20 )  » NUMB  OMC  20) «  T NRG IN ( 20  ) 
OIMENSIOf  MINBMBC20)  »  MINAPC20J 
COMMON /MARG IN/TMRG IN 
COMMON  PROSPGC12, 20,20) 

COMMON/F IXE  0/ RUNW 1 0 , OAHRA  0 , STA DEV, TOWID, PROCUT 
COMMON/GRP/ NSPG,NSPW,NGPW, STEP ,FRAI, ERUNWI.ERWID 
COMMON/E CUAL/MINAP,MINNBR,MINB MB, NUMBER 
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J10=0 

IFCtNGPW  .GT.12)  .AND.(PF0CUT.LT..99>>  J10*l 

IF ( (MGPH.GT. 12) .ANO. I PROC UT.LT . . 99) >  P  ROCUT *  PROCUT ♦•01 

EPSILN=fl. 

N=1 

A INLOC ( 1I=T0WID*0AMRAC 
Nl=RUNHIC/<  T0WI0*0AMRAD> 

X2=RUNWI0/(  TOWID+DAMRA  I) 

IF  I X2 • LT  .1. )  N=i 

IF (X2. LT  .1.  )  GO  TO  11 

IF  (1X2. LT. 2.)  .AN0.IX2.GE.1.))  N*1 

IFUX2.LT.2.)  .AND.  (X2.GE.1.)  )  GO  TO  11 

IFIX2.EQ.2.)  GO  TO  11 

IFIN1.EQ.2)  N  =  2 
IF(N1.EQ.2)G0  TO  11 
DO  10  1  =  2, N1 
X=I-1 

AIMLOCI I) =AIMLCC( 1)+X*(T0HI0»2.*0AHRAD) 

IF (AIMLOC(I) »GT.RUNWIO)AI MLOC(  I> ■ 

CAIMLOC  ( 1)  ♦X*|TOWIO*2.*OAMRAO)*TOHIO*OANRAD 
IF! AIHLOCII) .GT.  RUNW 10) A IMLOC ( I>  =0. 

IF  I AI MLOC (I ) • GT. RU KWID)GO  TO  11 

n«nh 

10  CONTINUE 

11  CONTINUE 
NUMAI H  =  N 
W=NUiA IH 

XiNUNAIM+1  •  --- 
Y* RUNWID/X 

- 2=W/2.  - - 

N8=H/2. 

Z1=Z-N8  . 

IF CZ1 . GT • 0* )  GO  TO  50 
DO  12  1*1, NO  ~ 

X1=I 

AIML3CI1  )=X1*Y 

12  CONTINUE 

N9*N8+ 1  ' 

DO  60  I=N9,N 
X1=M-I*1 

AIHLOC II) *RUNWI0-X1*Y 
60  CONTINUE 
GO  TO  70 

50  CONTINUE 
00  51  I*  1, N  8 
X1=I 

AIHLOC  ( I)  *X1*Y 

51  CONTINUE 
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N9=N8*2 
00  52  I=N9»  N 
Xi=N-I*t 

AIMLOC I I)=RUNWID-X1*Y 
52  CONTINUE 

CENTER=RUNWI0/2. 

AIMLOC  IN  8+1)  =C£NTER 
70  CONTINUE 

00  13  1=1, NUMAIN 
NUM30M  I  I)=l 

13  CONTINUE 
IFINSPW.GT.12)  SO  TO  20 
CALL  PRO (ALI NUMA  IM,  AIHLOC  ,  NUMB  OH) 

CALL  UNION!  NUMSOM,  NUMA  IM,  EPSIL  N,  TOTAL) 

GO  TO  21 

20  CONTINUE 

CALL  S INI NUMA IM, A INLOC , NUMBOM, TOTAL) 

21  CONTINUE 
MIN=NUMA IM 
PRONIN =1 -TOTAL 
IFIPROMIN.GE.  PROCUT)  MAX=  MIN 

- IFIPROMIN.GE. PROCUT)  GC  TO  100 

00  14  1=1,50 
00  15  J= 1 ,NUM AIM 
NUMBOM I J)=NUMBOM| J) >1 
IFINSPW.GT.12)  GO  TO  30 
CALL  P ROC ALI NUMA IM, A IMLOC , NUMBOM) 

CALL  UNION!  NUMBOM,  NUMAIN,  EPSILN,  TOTAL) 

GO  TO  31 

30  CONTINUE 

CALL  S  IMI  NUMA  IM,  A  IMLOC  , NUMBOM,  TOTAL) 

31  CONTINUE 
NUMBOM I J) sNUMBOMI J) -1 
TMRSINI J)=l. -TOTAL 

"15' "  CONTINUE 
N=  0 
X*0, 

00  16  K= 1, NUMAIN 
IFIX.LT .TMRGINIK) )  N=N*1 
1FIX.LT.TMRGINIK)  )  X=TMRG  INI  K) 

16  CONTINUE 

NUMBOM I N)=NUMBOM|N) *1 
IFITMRGIMN)  .GT. PROCUT)  GO  TO  100 
IFd.EQ.50)  PRINT*  *"  NOT  ENOUGH  TOTAL  WEAPONS" 
IF!  I.  E  Q.  50)  GO  TO  200 

14  CONTINUE 

100  CONTINUE  ' 
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N9=N8*2 
00  52  I=N9,N 
X1=N-I*1 

AIMLOC (I) =RU  NWID-X1*Y  - 

52  CONTINUE 

CENTER=RUNWI0/2. 

AIMLOC  !N  8*1 ) =  C£NT  ER 
70  CONTINUE 

00  13  1=1, NUMAIM 
NUM30M! I)=l 

13  CONTINUE 
IFCN3PW.GT.12)  SO  TO  20 
CALL  PRO  CAL< NUMA IM, AIHLOC , NUMB  ON) 

CALL  UNION!  NUMBOM, NUMA  IM,  EPSIL  N*  TOTAL) 

GO  TO  21 

20  CONTINUE 

CALL  SIM!  NUMAIM,  AIMLOC, NUMBOM,  TOTAL! 

21  CONTINUE 
MIN=NUMA IM 

PROMIM  =  1-T0T  AL  - 

IF (PROMIN.GE • PROCUT)  MAX=  MIN 

- - IFCPROMIN.GE.PROCUTI  GC  TO  100 

00  14  1=1.50 
00  15  J=1 .NUMAIM 
NUMBOM !  J ) =NUMBOM ( J) +1 
IFIN3PW.GT.12)  GO  TO  30 

.  CALL  PROCAL!  NUMAIM,  AIMLOC,  NUMBOM) 

CALL  UNION!  NUMBOM,  NUMAIM,  EPSIL  N,  TOTAL) 

GO  TO  31 

30  CONTINUE 

CALL  S  IM!  NUMA  IM,  AIMLOC  ,NUH80M,  TOT  AL) 

31  CONTINUE 
NUMBOM! J) =NUMBOM! J) -1 
TMR3IN ! J) =1 • -TOTAL 

"15'  “  CONTINUE 
N*  0 

“  x*o.  ~ 

DO  16  K=l, NUMAIM 

IFCX.LT  «TMRGIN!K)  )  N=N«-1  - 

IF!X.LT.TMRGIN!K)  )  X=T MRG  IN!  K) 

16  CONTINUE 

NUMBOM !N)=NUMBOM!N) ♦ 1 
IFUMRGIMN)  .GT.  PROCUT)  GO  TO  100 
IFCI.EQ.  50)  PRINT*,**  NOT  ENOUGH  TOTAL  WEAPONS 
IF! I*  E  Q*  50)  GO  TO  200 

14  CONTINUE 

100  CONTINUE  ‘  ' 
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MINN3  R=0 

00  10  3  I  sit  NUMAIM 
MINBM0 { I) =NUMB0M ( I) 

MINAP(  I)  =  AIHL0C<  I ) 

MINN3R=M INN8R*MINBMB< II 
103  CONTINUE 
N*0 

00  17  I=1,NUMAIM 
N=N*NUMB0M< I) 

17  CONTINUE 
MAX-N 

200  CONTINUE 

IF(JlO.EQ.l)  PROCUT=PROCUT-.01 

RETURN 

END 
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Subroutine  COUNT 


SUBR3UTINE  COUNT! ICOUNT,MINNUM» 


C 


THIS  SU3R3 UTINE  !COUNT)  ALLOWS  ONE  TO  VARY  THE  NUMBER 
OF  WEAPONS  TARGETEO  FOR  TH  f  AIM  POINTS.  IT  IS  USEO 
IN  SUBROUTINE  STRAT  TO  DETERMINE  THE  BEST  AIM  POINT 
STRATEGY  FOR  A  GIVEN  NUM8ER  OF  AIM  POINTS  WITH 
ICOUNT  WEAPONS  TARGETED  FOR  EACH  AIM  POINT,  I.E. 

EVERY  AIM  POINT  IN  STRAT  HAS  THE  SAME  NUMBER  OF 
WEAPONS  TARGETEO  FOR  IT. 

THE  FOLLOWING  VARIABLES  ARE* 

1.  MINNBR=*THE  CURRENT  UPPER  BOUND  ON  THE  NUMBER 
OF  WEAPONS  REQUIRED  TO  MAKE  THE  CUT  WITH  THE 
DESIRED  PROCUT 

2.  MIN  AP=AN  ARRAY  OF  THE  LOCATIONS  OF  THESE  AIM  POINTS 
WHICH  ACHE IV E  THE  CURRENT  UPPER  BOUNO  ON  THE 
NUMBER  OF  WEAPONS  REQUIRED 

. 3.  MINBMB=AN  ARRAY  OF  THE  NUMBER  OF  WEAPONS  TARGETEO 

FOR  EACH  AIM  POINT  IN  MINAP 

4Y  PROB=THE  CURRENT  PROBABILITY  OF  THE  -BEST-  AIM  POINT  STRATEGY 

5.  TMRGIN^AN  ARRAY  COMPUTED  IN  SUBROUTINE  BEST  THAT  HAS  THE 
CURRENT  PROBABILITIES  FOR  EACH  AIM  POINT  GROUPING 


C 


DIMENSION  MINAP!  20 1  ,  APHIN !  20 1  ,  HI NBMB! 20)  , MINBOM( 20) 
DIMENSION  TMRGTNI20)»V  ALUE !  20)  ,NUMBOM(20I 
DIMENSION  BESTAP!  20,20) 

DIMENSION  INDC20) 

DIMENSION  B VALUE! 20) 

COMMON / M  ARG IN/TMRG IN 

COMMON/OPT/BESTAP, INQ, BVALUE 

COMMON /EQUAL/ MINAP, MINNBR, MINBMB, NUMBER 


C 
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10  CONTINUE 

ICOUNT= ICOUNT  *1 
X* ICOUNT 
Y*MINNUN 
MAX=MINNUM 
ISTOP=Y/X 

IFIISTOP.LE.l)  SO  TO  100 

CALL  STRATt ICOUNT , MAX, VALUE,  NSET) 

CALL  BEST  Cl  COUNT  ,  MIN  NS  R,  MAX*  NUMB  ON,  NSET,  APMIN*  HINBOM,  MINNUM,  NGP1 

C  J45  IS  A  COUNTER  WHICH  SETS  THE  INITIAL  VALUE  OF  PROB, 

C  OEPENOINS  ON  WHETHER  ICOUNT  STARTS  AT  1  OR  2  IN  THIS 
C  SUBROUTINE 

IF(MINNBR.EQ. MINNUM)  60  TO  60 
J45=0 

IF<I30UNT.EQ.l)  J45=l 
IFCJ45.E0.1)  PROB=0. 

IF  (( ICOUNT.  EG.2) ■ AND. ( J45 •  EQ. 0)  )  PROB=0. 

IFCPROB.LT.TMRGIN(NGP)  »  60  TO  50 
60  TO  10 

50  CONTINUE 

PROB*  T  MRG IN  «  NGP) 

60  CONTINUE  ‘  .  .  . . 

NUMBER=NGP 

. MINNJM  =  0  *'■  “  . . . 

00-51  1=1,  NGP 

MINNUM  =  MINNUM«-MINBOM<I)  ""  . . .  - - -  - 

MINAP( I)=APMIN(I) 

MINBMB 1 1) =MINBOM( I)  “  “  - - 

51  CONTINUE 

IFCMINNUM.EQ.MINNBRI  GO  TO  IOO  - * - - 

GO  TO  10 

“100  CONTINUE  ~  .  '  ' 

RETURN 

Ei 
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Subroutine  STRAT 


SUBROUTINE  STRAT( ICOUNT,MAX,VALUE,NSET) 


THIS  SUBROUTINE  (STRAT  )  DETERMINES  THE  BEST  AIM  POINT 
STRATEGY  FOR  A  GIVEN  NUMBER  OF  AIM  POINTS  ANO  THE  FACT 
THAT  THESE  AIM  POINTS  HAVE  THE  SAME  NUMBER  OF  WEAPONS 
TARGETEO  FOR  EACH  OF  THEM.  THESE  STRATEGIES  ARE 
SAVED  IN  THE  BESTAP  MATRIX  ANO  THEN  ARE  PASSED  TO 
SUBROUTINE  BEST  FOR  DETERMINING  WHICH  OF  THE  THESE 
AIM  POINT  GROUPS  IS  ACTUALLY  THE  BETTER  AIM  POINT 
STRATEGY 

THIS  SUBROUTINE  WORKS  ON  THE  IDEA  THAT  IF  ALL  AIM 
POINTS  EXCEPT  THE  MIOCLE  TWO  ARE  HELD  FIXED  THE 
RESULTING  PROBABILITY  DISTRIBUTION  SHOULD  BE 
CONVEX.  THUS  WHEN  THE  PROBABILITY  OF  NOT  MAKING 
THE  CUT  BEGINS  TO  INCREASE  ONE  HAS  PASSED  THE  BEST 
LOCATION  FOR  THIS  CONFIGURATION  OF  AIM  POINTS 
WITH  AN  ODD  NUMBER  OF  AIM  POINTS,  THE  CENTER  AIM 
POINT  IS  ALSO  HELO  FIXED 

THE  FOLLOWING  VARIABLES  ARE  I  ~  . .  ~  - - 

1.  EPS IL N=A  CONVERGENCE  CRITERIA  FOR  SUBROUTINE  UNION 

WHICH  IS  SET  TO  ZERO  FOR  THE  USE  IN  THIS  THESIS  . 

2.  NSET=THE  CURRENT  NUMBER  OF  AIM  POINT  GROUPINGS 

"3.  CEN A IM=THE  CENTER  OF  THE  ACTUAL  RUNWAY  (CORRESPONDS  TO  THE 
HIOOLE  AIM  POINT  FOR  AN  ODO  AIM  POINT  GROUP) 

~5.~  AIMGP=AN  ARRAY  INDEXED  BY  AIM  POINT  GROUPING  AND  THE  AIM 
POINT  NUMBER  WITHIN  THAT  SPECIFIC  GROUPING 

6.  “  NSTEP=THE  NUMBER  OF  STEPS  (BY  FIVE  FEET)  FROM  THE  INITIALIZED 

POSITIONS  OF  THE  INNER  MOST  PAIR  OF  SYMMETRIC  AIM  POINTS 

7.  VAL  UE=AN  ARRAY  OF  THE  CURRENT  PROBABILITIES  OF  NOT  MAKING 
HE  CUT.  WHENEVER  THIS  PROBABILITY  DECREASES  FOR  A  SPECIFIC 
AIM  POINT  GROUP,  THE  NEW  PROBABILITY  IS  PLACED  IN  VALUE  AND 
HEN  THIS  PROBABILITY  IS  COMPARED  TO  THE  CURRENT  PROBABILITY 
STORED  IN  8VALUE.  IF  VALUE  IS  LESS  THAN  BVALUE,  THEN  THIS  AIM 
POINT  STRATEGY  IS  STORED  IN  BESTAP  BY  AIM  POINT  GROUP  NUMBER 
ANO  THE  LOCATION  OF  EACH  MEMBER  OF  THIS  GROUP 


76 


c 


DIMENSION  AIMLOCI20)*  NUMB OM( 20) «  AIMGPt  20*  20) » I fDX ( 20) » 
CBESTAP (20*20) « VALUE (20) ,INO( 20) t B VALUE ( 20) 

COMMON  PROSPGI12,20,20 
COMMON/A  IMG  RP/AIMGP 

COMMON/FIXE  D/RUNW I  D,OAPRAO»STA  DEV*  TOHIO,  PROCUT 
COMMON/GRP/NSPG*NSPW,NGPW,STEP,FRAI,ERUNWI*ERWIO 
COMMON/OPT/ BESTAP,  INO.EVALUE 


C  DETERMINE  THE  UPPER  BOUNO  ON  THE  NUMBER  OF  AIM  POINTS 

EPSILN=0. 

NSET=MAX/ICOUNT 

^  DETERMINE  THE  “BEST"  AIM  POINT  STRATEGY  FOR  EACH  AIM 
C  POINT  SET 

K*0  . 

10  CONTINUE 

K*K*1  '  ‘  -  - - - 

IF(K.GT.NSET)  GO  TO  1000 

' CENAIM=RUNM ID/2,  . 

IFfK.GT.  1)  GO  TO  20 

AIMGP(l,l)sCENAIM  '*  . . . . 

GO  TO  10 

“20*"  CONTINUE  “  -  - 

C  DETERMINE  WHETHER  ONE  IS  OEALING  WITH  AN  ODD  OR  AN  EVEN 
C  AIM  POINT  SET 

X*K 

N*X/2. 

Y*X/2 , 

Z*Y-N 

IFtZ.EQ.  Oi  GO  TO  30 
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Z  INITIALIZE  AN  000  AIM  POINT  SET 

AIMS?  <  K, N*1 ) =CENA IM 
N1*N*1 
00  21  1*1.  N 
X=  I 

AIMG^IK.  I)=X*0AMRA0 

21  CONTINUE 
M1*N*2 

00  22  I=M1.  K 
X* 

AIMGPIK.  1>  =  RUNWI0-X*DAMRAD 

22  CONTINUE 
GO  TO  40 

30  CONTINUE 

C  INITIALIZE  A  A  EVEN  AIM  POINT  SET 


N1=N 

00  31  1*1.  N 
X*I 

AIMGPIK. I)=X*OAHRAO 

31  CONTINUE 
M1*N*1 

00  32  I*H.K 
X«K-I*1« 

AIMGP(  K,  I)  =  RUNHIO-X*OAMRAO 

32  CONTINUE 
40  CONTINUE 

THE  ITERATIVE  SHEME  TO  DETERMINE  THE  “BEST**  AIM  POINT 
STRATEGY  PER  AIM  POINT  SET 

- DO  41  1=  1.  K 

AIMLOC 1 1 ) =A IMGPI K  » II 
- -  NUMBON(I)  =ICOUNT 

"41  CONTINUE 
NUMAIM=K 

IFINGPW.GT.12)  GO  TO  50 

CALL  PRO CAL  I NUMA IN »  AIMLOC . NUMB0M1 

CALL  UNION!  NUMBOM.  NUMAIM.  EPSILN,  TOTAL) 

GO  TO  51 

50  CONTINUE 

CALL  SIMINUMAIM, AIMLOC.NUMBOM,  TOTAL) 

51  CONTINUE 
VALUE I  K)  =TOTAL 
00  42  Is  l.K 

BESTAP  IK  » I)  =  A IMGP IK. I ) 

42  CONTINUE 
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J36=<- 1 
00  SO  IS1«J36 

OIS=AIMGP«, I*1)-AIMGP<K» I) 

IFTOIS.LE.OAMRAO)  GO  TO  70 
60  CONTINUE 
GO  TO  60 

70  CONTINUE 

00  71  Is ly K 
BESTAP  IK*  I)  =  0  . 

71  CONTINUE 
GO  TO  10 

80  CONTINUE 

IFIN1.EQ.N)  GO  TO  200 

C  WORKING  WITH  THE  000  AIM  POINT  SETS 

L*N1-1 

CL  IS  AN  INDEX  FOR  DETERMINING  WHICH  PAIR  OF  SYMMETRIC 
C  AIM  POINTS  IS  BEING  MOVED*  THE  SMALLER  L  IS,  THE  CLOSER 
C  ONE  IS  TO  THE  MIOOLE  PAIR«  WITH  L* 1  BEING  THE  MIDDLE  PAIR. 

00  101  1*1,  L 
INO<I>=I 
101  CONTINUE 

BVALUE  <K) =1 • 

Ll=l 

'*105  CONTINUE 

IF(Ll.GT.L)  GO  TO  10 
110  CONTINUE 
J21=0 


C  J21  IS  A  COUNTER  ON  THE  NUMBER  OF  TIMES  THE  AIM  POINT 
"C  STRATEGY  EXCEEOS  THE  DIFFERENCE  CRITERIA 

•*  "  "  T*L1  .  . - . 

VALUE  00*1. 

-  IF(K.eq.  :r~co  to  i5o 

AIMGPf  K, N1-IN0<I>  »  *AIMGP(  K*N1-IN0CI> )*5« 

AIM3PC  K,  NOINO(I))  =  AI  MGPt  K,  Nl*  INO(  I>  )  -5. 

C  MOVES  EACH  OF  THE  AIM  POINTS  DESIGNATED  BY  LI  <I»  ^  FEET 
C  INWARD 
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Jl*Nl-INOCI)  *1 
J2*Nl*  1 
J4=N1-1 


C  INITIALIZES  ALL  AIN  POINTS  INBOARD  OF  THE  TWO  JUST  MOVED 
C  (ALL  PAIRS  WITH  L1<I> 

00  112  M2=J1,J4 
X*X*1. 

AIMGPC  K»  M2)  -AIHGPC  K,N1-IND(D)  *X  *OAMRAO 

112  CONTINUE 

IFCAIMGPCK,  J4>  .GE.CENAIM)  Ll*Ll*i 

IFCAINGPCK,  J4)  .GE.CENAIM)  GO  TO  105 

IF  ((GENA  IN- A  IMGP(  K,  J4I  I  .LT.QANRAO)  L1=L1*1 

IFCCCENAIM-AIMGPCK,  J4)  >. LT.QANRAO)  GO  TO  105 

J3=N1*INC(I)  -1 

DO  113  N2=J2,J3 

X*K-M2 

AIHGPC  K,  M2)  = AIMGPC  K,  IN  TCI)  *N1)  -X*OAMRAO 

113  CONTINUE 
GO  TO  150 

120  CONTINUE 
1*1*1 

IF  Cl. GT, L)  GO  TO  10 
GO  TO  110 

150  CONTINUE 

NSTEP=(CENAIN-AINGP(K,Nl-l))/5. 

1*1 

151  CONTINUE 

1*1*1  . 

IF  Cl.  EQ.  2)  GO  TO  155 
IFCI.GT.NSTEP)  GO  TO  154 
AIMGPC  K,N1-1)*AIMGP(K,N1-1>  *5. 

C  CHECK  TO  MAKE  SURE  THAT  THE  CENTER  PAIR  OF  AIM  POINTS 
C  Cll  HAVE  NOT  CROSSEO  THE  CENTER  POINT  OF  THE  RUNWAY 
C  C2I  ARE  NOT  LESS  THAN  ONE  DAMAGE  RAOIUS  FROM  THE  CENTER 
Z  OF  THE  RUNWAY 

IF CAIMGP ( K,  N 1-1)  .GT.CENAIM)  Ll*Ll*l 
IFCAIMGP  CK,N1-1>  .GT.CE)AIM)  GO  TO  105 
IFCCCENA  IM-AIMGPC  K,N1-1))  .LT .0 AMRAO)  Ll*Ll*l 
IFICGENA  IH-AIMGPCK,Nl-in  .LT.OAMRAD)  GO  TO  105 
AIMGPC  K,Nl*l)=AIMGTCK,Nl*ll-5. 

155  CONTINUE 

00  152  11*1,  K 
AIMLOCC  Ii»  =  AIMGPCK,Il» 

152  CONTINUE 
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C  ALLOWS  THE  FIRST  OCCURRENCE  OF  OIF  EXCEEDING  IT'S  CRITERIA 
C  TO  NOT  STOP  THE  ITERATIVE  PROCESS.  THE  SECONO  TINE  OIF  EXCEEOS 
C  ITS  CRITERIA,  THE  CURRENT  fiASE  PAIR  OF  AIN  POINTS  (LI)  IS  MOVE 
C  5  FEET  INBOA  AO  ANO  THE  ITERATIVE  PROCESS  STARTS  AGAIN 

IF(I.  GT.  2)  CALL  APROX C  AIHLOC,  NUMAIM,  NUM80N*  OIF ) 

IF(I.LE.3)  niF=.m 

IFiOIF.LT.O.)  PRINT*, "STARTED  TOO  CLOSE  TO  THE  EDGE" 

IFIDIF.G  7.. 05)  J21=J21*1 

IF((J21.GT.l)  .AND. (K.EQ.3))  L1=L1*1 

IF(( J21.GT. 1) .AND. ( K.EQ.3) )  GO  TO  105 

IF( J21 .GT .1)  GO  TO  110 

IF (NGPW. GT .12)  SO  TO  160 

CALL  PROCALI NUMA IN, AIHLOC , NUH80M) 

CALL  UNI GN(  NUMBOM, NUMA IN,  EPSILN, TOTAL! 

GO  TO  161 

160  CONTINUE 

CALL  S IN ( NUMA IN, A IMLOC , NUMBOM, TOTAL) 

161  CONTINUE 

IF((NGPW.GT.12)  .ANO.(I.EQ.2)  )  GO  TO  163 

CHECK  TO  SEE  IF  TOTAL,  THE  PROBABILITY  OF  NOT  MAKING 
A  CUT,  IS  STILL  DECREASING.  IF  IT  IS,  RESET  VALUE  (K)  TO 
TOTAL.  IF  IT  IS  NOT,  THEN  MUST  MOVE  THE  CURRENT  PAIR  (LI) 

5  MORE  FEET  INBOARO  ANO  ITERATE  AGAIN. 

IF((VALUE(K)  .LT. TOTAL)  .ANO.  ( K. ED  .3)  )  L1=L1»1 
IF((VALUE(K)  .LT. TOTAL)  .AND. (K.EQ.3))  GOTO  105 
IF((VALUE(K) .LT. TOTAL) .AN O. C I. ED .31 )  L1=L1*1 
IFUVALUE(K)  .LT. TOTAL)  .AN  O.  ( I.  ED  .  3)  )  GOTO  105 
IF(VALUE(K)  .LT. TOTAL)  GO  TO  110 

163  ’  CONTINUE  -  -  ~  . 

VALUE ( K)  =TOTAL 

IF (VALUE (K) .GT.BV ALUE( K) )  GO  TO  151 
BVALUE (K) =V ALUE(K) 

'  DO  153  J10=l,  K  .  “  " 

BESTAP (K, J10!=AIMGP(K, J10) 

‘  153  CONTINUE 
GO  TO  151 
156  CONTtNUE 

PRINT*,-  IS  CONTINUALLY  DECREASING  ON  THIS  INTERVAL.  K*  -,K 
L1=L1*1 
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200  CONTINUE 

C  WORKING  WITH  THE  EVEN  AIH  POINT  SETS  - 

L*N1 

L  IS  AN  INDEX  FOR  DETERMINING  WHICH  PAIR  OF  SYMMETRIC 
AIM  POINTS  IS  BEING  MOVED.  THE  SMALLER  L  IS*  THE  CLOSER 
ONE  IS  TO  THE  MIDDLE  PAIR,  WITH  L=  1  BEING  THE  MIOOLE  PAIR. 

00  201  1=1*  L 
INOUI  *1 

201  CONTINUE 
BVALJEt  K»*l. 

Ll*l 

205  CONTINUE 

IFCL1.GT.L)  GO  TO  10 
210  CONTINUE 
I*L1 

VALUE  <Kl=l. 

C  J22  IS  A  COUNTER  ON  THE  NUMBER  OF  TIMES  THE  AIM  POINT 
C  STRATEGY  EXCEEDS  THE  DIFFERENCE  CRITERIA 

J22=0 

IFIK.EQ.  21  GO  TO  250 
IF  (LI.  ED.  1)  GO  TO  250 

C  MOVES  EACH  OF  THE  AIH  POINTS  DESIGNATED  BY  LI  Cl»  5  FEET 
C  INWARD 

AlM5PCK,Nl»l-INOIinaAIMGPCK,Nl»l-INDf  I))»5. 

AIMGPI  K,N1*  INOm  )  =AIMGPC  K,  Nl»  IND<  III  -5. 

*"  J1=N142-IN0<  II 
IF(Ll.EQ.l)  J1=N1 
X*Q. 

C  INITIALIZES  ALL  AIM  POINTS  INBOARD  OF  THE  TWO  JUST  MOVEO 
C  tALL  PAIRS  WITH  L1<I> 

00  212  M2=J1,N1 
X*X*1. 

AIMGP t  K*  M2I =AIMGP(K*N1*1-IND(III  ♦X*OAMRAO 

212  CONTINUE 

IF(  AI MGP <K»  Nil  .GE.CENAIM)  L1=L1*1 
XF(AIMGP(K,N1I .GE.CENAIMI  GO  TO  205 
J2*Nt*l 

J3«N1*IN0<I> -1 
IF  (LI.  EQ.  1)  J3=N1  ♦  1 
00  213  M2aJ2*  J3 
X*K-M2 

AIMGPf  K* P2) *AIMGP(K* IND<I> ♦Nil -X •OAHRAO 

213  CONTINUE 


82 


onn 


CHECK  TO  MAKE  SURE  THAT  TH  t  CENTER  PAIR  OF  AIM  POINTS 
< 1 )  HAVE  NOT  CROSSED  THE  CENTER  POINT  OF  THE  RUNWAY 
121  ARE  NOT  LESS  THAN  ONE  DAMAGE  RADIUS  APART 

IF  IAIMGP  IK,  N 1)  .GE.AIMGFIK,  J2II  Ll=Ll*l 
IF! AIMGP  IK,  Nil.GE.  A  IMG  P(K  «  J2) I  SO  TO  205 
OIS=AIMGPIK,  J2)-AIMGPIK,Ni) 

IF(DIS  .L  E.OAMRADJ  L1=L1*1 
IFIOIS.LE.OAMRAO)  GO  TO  205 
GO  TO  250 
220  CONTINUE 
1=1*1 

IFII.GT.L)  GO  TO  10 
GO  TO  210 

250  CONTINUE 

NSTEP* I  CENA IM-AIMGP(K,Nll  1/5. 

1=1 

251  CONTINUE 
1=1*1 

IFII.EQ.  2)  GO  TO  255 
IFCI.GT.NSTEPI  GO  TO  254 
AIMG®  I  K,  N1)=AIMGP(  K,N1)*5. 

AIMGPI K,N1*1»  =  AIMGP<K,M*l>-5. 

255  CONTINUE 

J=N1* 1 

_ IF<AIMGP(K,  Nil  .GE.AIMGP«K,N1*1»)  L1=L1*1 

IF(AIHGPIK,N1)  .GE.AIMGFIK, Nl*ll)  GO  TO  205  . .  ' 

__  0IS=AIMGPIK,N1*1) -AIMGFlKtNl) 

IFIOIS.LE.OAMRAO)  Ll=H*i  - - -  -  - - - 

IFCDIS.LE.DAMRAOI  GO  TC  205 

. OO  252  J*l,  K  -  ... 

AIMLOC ( J) =AIMGP<K, J> 

"252  CONTINUE  - -  -  -  - - 

C  ALLOWS  THE  FIRST  OCCURRENCE  OF  OIF  EXCEEDING  IT'S  CRITERIA 
C  TO  NOT  STOP  THE  ITERATIVE  PROCESS.  THE  SECOND  TIME  OIF  EXCEEDS 
C  ITS  CRITERIA,  THE  CURRENT  BASE  PAIR  OF  AIM  POINTS  ILi)  IS  MOVE 
C  5  FEET  IN30AR0  AND  THE  ITERATIVE  PROCESS  STARTS  AGAIN 

IFd.ST.3l  CALL  APROX I  AIM LOC, NUM  AIM,  NUHBOM,  DIF) 

IFII.LE.3I  0 IF=.  01 

IFIOIF.LT. 0.)  PRINT*, -STARTED  TOO  CLOSE  TO  THE  EDGE* 

IFCDIF.GT..05)  J22=J22*l 

IF (I J22.GT.  1)  .AND.  (K.E0.2I)  Ll*Li*l 

IFII  J22.GT.  1) .AND.  IK.EQ.2I)  GO  TO  205 

IFIJ22.GT.il  GO  TO  210 

*  • 
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IF<NSPW.GT.121  GO  TO  260  ' 

CALL  P  ROCAL I NUMA IM  ,A IMLOC *  NUMBOM) 

CALL  UNICNlNUM60M,NliMAIHt  EPSILN,  TOTAL) 

GO  TO  261 

260  CONTINUE 

CALL  S  IMCNUMAIM,  A  IMLOC  .NUMBOM,  TOTAL! 

261  CONTINUE 

IFUNGPM.GT.12)  .  ANO.  ( I  .Ed  .  2!  !  GO  TO  263 

C  CHECK  TO  SEE  IF  TOTAL*  THE  PROBABILITY  OF  HOT  HAKING  A  CUT. 
C  IS  STILL  DECREASING.  IF  IT  IS,  RESET  VALUECK!  TO  TOTAL. 

C  IF  IT  IS  NOT,  THEN  MUST  MOVE  THE  CURRENT  PAIR  (LI!  5 
C  MORE  FEET  INBOARD  AND  ITERATE  AGAIN. 


IF <(VALUECKI.LT. TOTAL!  .AND. <K. EQ  .21 1  L1=L1+1 
IF(<VALUE<K)  .LT.TOTAL1  .AND  .  <  K.  EQ  .2!  !  GO  TO  205 
IF<<VALUE<K1  .LT. TOTAL!  .AND. CI.EQ. 3!!  L1=L1  +  1 
IF<<VALUE<K1  .LT.  TOTAL!  .ANO. (I. EQ. 31!  GO  TO  205 
IFCVALUECK! .LT. TOTAL!  GO  TO  210 
263  CONTINUE 

VALUE <  K1 =TO T AL  -  - 

IFCVALUECK!  .GT.BV  ALUE  <  K!  )  GO  TO  251 

BV ALUE ( K1 =V  ALUE<  K1 

OO  253  J  10=1,  K 

BESTAP (K»J101  =  AIMGP(K, J10 1 

253  CONTINUE 

-  60  TO  251  '  '  .  . . 

254  CONTINUE 

- PRINT*,"  IS  CONTINLALLY  DECREASING  ON  THIS  INTERVAL 

L1=L1*1 


GO  TO  205 
1000  CONTINUE 

- RETURN 

END 


•  K* 
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Subroutine  APROX 


SUBROUTINE  APROX( AIMLOC.NUMAIH,  NUMBOM,DIF) 


C 


THIS  SUBROUTINE  (APROX)  COMPARES  THE  PROBABILITY  OF  DAMAGING 
THE  MIDDLE  GROUP  WITH  THE  PROBABILITY  OF  OAMAGING  THE 
FIRST  GROUP.  THIS  ALLOWS  FOR  A  QUICK  CHECK  ON  THE  AIM 
POINT  STRATEGY  UNOER  CONSIDERATION  IN  SUBROUTINE  STRAT. 

THE  FOLLOWING  VARIABLES  ARE  I 

1.  OAMLOC=AN  ARRAY  OF  THE  AIM  POINT  LOCATIONS  RELATIVE  TO  THE 
EFFECTIVE  RUNWAY  WICTH 

2.  CENTER=C ENTER  OF  THf  ACTUAL  RUNWAY  WIDTH 
'3.  5RPLE N=L ENGTH  OF  A  GROUP 

•».  HAFLEN=HALF  OF  GRPL  fN 

5.  MIOGRP=THE  STARTING  POINT  OF  THE  "MIDDLE*  GROUP 

6.  EQUPRO=AN  ARRAY  OF  THE  SECTION  BOUNDARIES  IN  TERMS  OF 
STANDARD  DEVIATION  FROM  ITS  AIM  POINT.  THE  ARRAY  IS  INDEXED 

BY  GROUP  NUMBER  i WHERE  1  IS  THE  FIRST  GROUP  AND  2  IS  THE  MIDDLE 
GROUP),  AIM  POINT  NUMBER,  AND  SECTION  NUMBER  WITHIN  EACH  GROUP 

7.  MONOR= AN  IMSL  SUBROUTINE  THAT  RETURNS  PCX<A)  FOR  A  N(0,1> 

-  DISTRIBUTI0N 

S.  PROEQU=AN  ARRAY  OF  THE  PROBABILITIES  THAT  A  WEAPON  IMPACTS 
WITHIN  A  SECTION.  THE  ARRAY  IS  INDEXED  THE  SAME  AS  EQUFRO. 


e 


"  DIMENSION  E QUPRO(  2,  20, 25)  ,  S1I25) ,  PROCC  2)  , PR0EQUC2, 20,25) 
DIMENSION  DAMLOC(20), NUMB ONI  20  > , AIMLOCC20) 

. COMMON /F IXED/RUNW 10 , DA  7RA  D, STA  DEV, TO HID, PROCUT 

COMMON/GRP/NSPG, NSPW, NGPW , STEP, FRA  I, ERUNWI, ERWID 
COMMON/APPROX/EQUPRO,Sl,PROD,PROEQU 


AD-A094  807  AIR  FORCE  INST  OF  TECH  WRI6HT-PATTERSON  AFB  OH  SCHOO-- ETC  F/6  15/7 

A  GENERALIZED  COMPUTER  MODEL  FOR  THE  TARGETING  OF  CONVENTIONAL  — ETC(UI 
DEC  80  J  C  PEMBERTON 

UNCLASSIFIED  AFIT/60R/OS/80D-6  NL 


DO  1  1  =  1, NUMAIM  -  "  “ 

OAMLOC  ( I)  =AIMLOC(  II+DAFRAO 
CONTINUE 

CENTER  =RUNW ID/2. 

GRPLEN  =  F  RAI*STEP+T  CHID+OA  MRAO 
HAFL*N=GRPLEN/2. 

MIOGRP=CENTER-HAFL£N 
NSPS*1=NSPG  +  1 
00  2  1=1,2 
00  2  K  =  l, NUMAIM 
00  2  J=1,NSPGP1 
X=0. 

IF(I.EQ.  Z)  X=HIDGRP 

y=j-i 

IFII.EQ.NSPG)  Y=Y ♦FRAI 

EQUPRO (I, K, J) =<(Y*STEP) *X -OAMLOC <  K) > /STAOEV 

CONTINUE 

00  3  1=1,2 

00  3  K= 1, NUMAIM 

DO  3  J  =  1 ,  NSPG 

CALL  HONORCEQUPROCI,K,J>  , BONDL) 

L=Jd 

“CALL  MONORC  EQUPRO  C  I,  K,L),B0N0R1 

PROEQU  C  I  «K«  J)  =BONOR-BONDL 

CONTINUE 

00  6  1=1,2 

PR0D( I ) =1 . 

00  4  J*l,  NUMAIM 
S1CJ>  =0. 

00  5  K=1,NSPG 

SIC  J)  =S1IJ)  ♦PROEQU  Cl,  J,  K) 

CONTINUE  . 

S1CJI  =S  1 C  J)  **NUMBOMCJ> 

PROOC  11  =  FROOC  I)  *S1  ( J)  . . 

CONTINUE 

CONTINUE  "  . 

0IF=®R00C2) -PROOC  1) 

RETURN 

ENO 
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Subroutine  BEST 


SUBROUTINE  8EST< ICOUNT  ,MIN,MAX, NUMBOM, NSET, APMIN, MIN80M,  MINNUM, 
CNGPI 


C 


THIS  SUBROUTINE  (BEST)  TAKES  THE  MOST  "EFFICIENT"  AIM  POINTS 
STRATEGT  PER  AIM  POINT  GROUP  IBESTAPI  ANO  THEN,  THROUGH 
MARGINAL  ANALYSIS,  DETERMINES  WHICH  OF  THESE  AIM  POINT 
STRATEGIES  REQUIRES  THE  MINIMUM  NUMBER  OF  WEAPONS. 

THE  FOLLOWING  VARIABLES  A  PE* 

1.  APM IN= AN  ARRAY  OF  THE  LOCATIONS  OF  THE  CURRENT  "BEST"  AIM 
POINT  STRATEGIES 

2.  MINBON=AN  ARRAY  OF  THE  NUMBER  OF  WEAPONS  TARGETED  FOR 
EACH  SPECIFIC  AIM  POINT  IS  APMIN 

3.  EPS ILN=A  CONVERGENCE  CRITERIA  FOR  SUBROUTINE  UNION 
WHICH  IS  SET  TO  ZERO  FOR  THE  USE  IN  THIS  THESIS 

4.  NMA  X=THE  LATEST  UPPER  BOUND  ON  THE  TOTAL  NUMBER  OF  WEAPONS 

5.  TMRGIN=AN  ARRAY  USED  TO  COMPARED  THE  MARGINAL  INCREASE 

IN  THE  PROBABILITIES  OF  MAKING  A  CUT  GIVEN  THAT  ONE  WEAPON  IS 
ADDEO  TO  AN  AIM  POINT 

"6."  MINNUM=THE  TOTAL  NUMBER  OF  WEAPONS  IN  MINBOM 
7.  NSP=N6=C0UNTERS  TO  INDICATE  THE  "BEST"  AIM  POINT  GROUP. 


C 


DIMENSION  8ESTAP<  20,  20  )»NUM8  0M(  20)*APMIN<20>  ,MINBOMC20)  , 
CTMRGIN  ( 2  0)  ,  A  IMLOC  <  20) 

DIMENSION  IN0C2Q) 

DIMENSION  B V ALUE( 20) 

COMMON  PROSPG<12,20,20) 

COMMON/MARG  IN/TMRG IN 

COMMON/ FIXE  D/RUNW  IO, OAMRAO,STA DEV, TOWIO,PROCUT 
COMMON/GRP/ NSPG»NSPW, NCPW , STEP, FRAI, ERUNWI, ERWID 
COMMON/OPT/BESTAP,  INO,  6VALUE 


C 
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00  1  1*1,20 

TMRGIM  (11=0. 

APMIN (  I)  *0. 

MIN30M ( I ) =0 

I  CONTINUE 
EPSILN=0 • 

J21=0 

C  SETS  THE  CENTER  OF  THE  RUNWAY  AS  THE  AIM  POINT  IN  AIM  POINT  GROUP  1 

BEST  A  P  (1,1)  *RUNWI 0/2  • 

N6*0 

L*0 

C  L  IS  A  COUNTER  ON  THE  NUMBER  OF  THE  AIM  POINT  GROUP 
C  AFTER  AIN  POINT  GROUP  1  IS  CHECKED  ONCE  IT  NEEO  NOT  BE  RECHECKEO 

IFClCOUNT.Eq.il  TMRGIN(l)  =0. 

IFCICOUNT.Eq.il  GO  TO  10 
IFIICOUNT.EQ  .21  L*1 
NUHAIH=L 

IFCIGOUN  T.EQ.2I  GO  TO  20 
IFCICOUN  T.GT.2)  1*1 
10  CONTINUE 
L*L*1 
NUMAIH*L 

. . .  IFCL.GT.NSET)  GO  TO  110 

_ _ J8=l  .  . . .  . . . 

Z  SETS  THE  "BEST-  AIM  POINTS  FOR  EACH  AIM  POINT  GROUP, 

C  AS  DETERMINED  TN  STRAT  ,  INTO  AIMLOC  FOR  USE  IN  EITHER  PROCAL  OR  SIN 
C  SETS  THE  INITIAL  ALLOCATION  OF  WEAPONS  TO  EACH  AIM  POINT  AT  ICOUNT 

IFCBESTAFCL, 1) .EQ.O.I  GO  TO  10 

00  11  I*  J8,  L 

AIMLOC 1 1 ) =BESTAP ( L, 1 1 

NUMBOM (II =ICOUNT 

II  CONTINUE 

SETS  THE  MAX  NUMBER  OF  TOTAL  WEAPONS  ALLOCATED  TO  MAX-ICOUNT*L 
STOPS  WHEN  NO  MORE  WEAPONS  CAN  BE  ALLOCATED  WITHOUT  EXCEE0IN5  MAX 

‘  I*MAX-CL*ICOUNT) 

IFCI.EQ.  0)  GO  TO  30 


i 


88 


C  DETERMINES  THROUGH  MARGINAL  ANALYSIS  WHICH  AIM  POINT  LOCATION  WITHIN 
C  AN  AIM  POINT  GROUP  WILL  OECREASE  TOTAL  THE  QUICKEST,  GIVEN  THAT  ONLY 
C  ONE  MORE  WEAPON  IS  ADOEO  TO  ANY  OF  THE  AIM  POINT  LOCATIONS 

00  12  J=  it  I 

00  13  J1=J8,L 

NUM80M ( J i)  =  NUMBOM ( J 1 ) ♦  1 

IFCN3PW.GT.12)  GO  TO  50 

CALL  PROCALC NUHAIM, AIM LOC,NUMBOHI 

CALL  UNIONC  NUMBOH, NUHAIM,  EPSILN,  TOTAL) 

GO  TO  51 

50  CONTINUE' 

CALL  SIHC NUHAIM, A IHLOC, NUHBOH, TOTAL) 

51  CONTINUE 

NUMBOH (Ji)sNUMBOHC Ji)-i 
TMRGINCJ  1)  =  1. -TOTAL 

IFCNGPW.GT.12!  TMRGINC  Jl) *1.01-TOTAL 

13  CONTINUE 
N=0 

X=0. . .  '•  - 

00  14  K=J8,  L 

IFCX.LT.TMRGINCK) )  N=NM 
IF  CX.  LT«  TMRGINCK)  )  X=THRG  INC K) 

14  CONTINUE  '  ~ 

NUMBOH  C  N) =NUMBOMC  N) +1 

-  IF  CTHRGINCN)  «GE«PROCUT)  GO  TO  15 
12  CONTINUE 
GO  TO  10 

30  CONTINUE 

- - IF<NGPW.GT.12)  GO  TO  31  - *  - — . . . 

CALL  PROCALC NUMAIM,AIMLOC, NUMBOH) 

- - CALL  UNIONC  NUMBOH, NUMAIM,  EPSILN,  TOTAL)  - 

‘ ‘ GO  TO  32 

31  CONTINUE 

CALL  SIMCNUMAIM,AIMLOC,NUMBOM, TOTAL) 

32  CONTINUE 
TMRSINCLM1. -TOTAL 

IFCNGPW.  GT.  12)  TMRGINCL)  =  1.01-TOTAL 
IFCTMRGINCL) ,GE«PROCUT)  GO  TO  15 
GO  TO  10 

15  CONTINUE 

C  COUNTS  UP  THE  NUMBER  OF  TOTAL  WEAPONS  USEO 
NMAX*0 

00  16  K*  J  8,  L 
NMAX=NMA >*NUHBOHC K) 

16  CONTINUE 

IFCCTHRGINCL) .LT.TMRGINC1)).AN0.  CNHAX.EQ.MAX) .AN0.J21.EQ.il 
C  GO  TO  10 
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C  IF  THE  TOTAL  NUMBER  OF  WEAPONS  USEO  IS  LESS  THAN  THE  CURRENT  UPPER 
C  90UN0  THIS  AIM  POINT  STRATEGY  IS  STORED  IN  APMIN  ANO  MINBOM 
C  ANO  MAX  IS  NOW  SET  TO  THE  NEW  UPPER  BOUND*  NMAX 

IFCNMAX.LT.  MAX)  N6*L 
IFINMAX.LT. MAX)  MAX=NMAX 
IF  I  MAX .E Q.HIN)  GO  TO  100 
IF(NMAX.EQ.MAX)  N6=L 
IF  (NMAX.  EQ.  MAX)  SO  TO  100 
GO  TO  10 
20  CONTINUE 
J2 1=1 

NUM30M ( 1) =MAX 
AIML0GI1)=BESTAPC1,1) 

NUMAIMsl 

IF CNSPW. GT. 12)  GO  TO  70 

CALL  PROCALC  NUMAIM, AIMLOC  *  NUMBOM) 

CALL  UNIONC  NUMBOM*  NUMAIM*  EPSIL  N*  TOTAL) 

GO  TO  71 

70  CONTINUE 

CALL  S IMC  NUMA IM*  AIMLOC,  NUMBOM,  TOTAL) 

71  CONTINUE 
TMRGIN<1)=1. -TOTAL 

100  CONTINUE 

MINNJMrO 

J8*l 

OO  101  I=J8«L 
APMINC  I)  *BESTAP(L,I) 

MINBOM ( I) =NUMBOM(  I) 

MINNUM=MINNUM*MINBOMCI> 

101  CONTINUE 

..  .  ....  I(r ( NM A X •  E Q.  M I N)  GO  TO  110  ~ 

GO  TO  10 
110  CONTINUE 
NGP=N6 

— .  IFCNS.EQ.  0)  NGP*1  —  ‘  . . .  . .  . 

RETURN 

ENO 
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Subroutine  PROCAL 


SUBROUTINE  PROCAL ( NUM AIM, AIMLOC,  NUMBOM) 


THIS  SUBROUTINE  (PROCAL)  CALCULATES  THE  PROBABILITIES 
OF  A  WEAPON  LANDING  IN  A  SFECIFIC  SECTION  GIVEN  AN 
AIM  POINT  STRATEGY  (THE  NUMBER  OF  AIM  POINTS,  THEIR 
LOCATION,  AND  THE  NUMBER  OF  WEAPONS  TARGETED  FOR  EACH 
SPECIFIC  AIM  POINT).  THIS  INFORMATION  IS  STORED  IN 
PROSPG  BY  GROUP  NUMBER,  AI)  POINT  NUMBER,  AND  SECTION 
NUMBER  WITHIN  THIS  GROUP. 

THE  FOLLOWING  VARIABLES  ARE  * 

1 •  SECBO  N=AN  ARRAY  THAT  HAS  THE  BOUNDARIES  UN  TERMS  OF 
STANDARD  DEVIATIONS)  OF  THE  SECTIONS 

2.  SECFRA=AN  ARRAY  THAT  HAS  THE  BOUNDARIES  OF  THE  LAST 
SECTION  IN  A  GROUPING.  THESE  BOUNDARIES  ARE  NOT 
A  MULTIPLE  OF  STEP  AND  MUST  BE  COMPUTED  SEPARATEV 

3.  PROSEC=AN  ARRAY  WHICH  HAS  THE  PROBABILITIES  OF 
A  WEAPON  LANDING  IN  THAT  SECTION.  IT  USES 
AN  IMSL  SUBROUTINE  CALLEO  HONOR  TO  COMPUTE  THIS 
PROBABILITY. 

4.  PROFAC=AN  ARRAY  OF  THE  PROBABILITIES  A  WEAPON  LANOS  IN  THE 
SECTION  OF  A  SPECIFIC  GROUP.  IT  ALSO  USES  HONOR  (SEE  3) 


"  DIMENSION  A  I MLOC  (  2  0 )  ,  NUMBOM  (  20 )  ,  SECBON  (20,  20),  PROSEC  (  20,  20) 
DIMENSION  PROFAC(20, 20), OA  MLOC  (20), SECFRA  (20,20) 

COMMON  PROSPG(12, 20,20) 

COMMON/ JOHN/  SECBON, PROSEC, PROFA C, DA MLOC, SECFRA 
COMMON/FIXE O/RUNW ID, D AMRA D,STA OEV , TOWID, PROCUT 

common/grp/nspg,nspw,ngpw,step,frai,erunwi,erwid 
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00  1  1*1 ,NUMAIM 

0AHL0C  ( I)  =A  INLOCC I)  ♦OAPRAO 

00  1  J=1,NSPW 

X*J-1 

SECBON  Cl,  J)  *  ( X*STEP-0A  M.0CCIH  /STAOEV 
1  CONTINUE 

DETERMINE  THE  PROBABILITY  OF  A  WEAPON  IMPACTING  WITHIN 
EACH  SECTION  FOR  EACH  AIM  POINT 

00  2  1*1 ,NU HA IM 
NSPWM1=NSPW-1 
00  2  J  =  l, NSPHM1 

NEED  TO  COMPUTE  PROBABILITY  FROM  -INF  TO  LEFT  BONO  CBONDL) 
CALL  HON  CRC  SECBON ( I, J) , BO NOLI 

NEEO  TO  COMPUTE  PROBABILITY  FROM  -INF  TO  RIGHT  BONO  CBONDR) 
K*  J*1 

CALL  HON  CRC  SECBON Cl, K) ,80 NOR) 

PROBABILITY  FOR  THE  SECTION  IS  THE  DIFFERENCE  BETWEEN  THE 
THE  LEFT  BONO  AND  THE  RIGHT  BOND  CPROSECI 

PROSEC ( I , J) =BONOR-BONOL 
CONTINUE 

CALCULATE  THE  FRACTION  OF  THE  TOWIO  NOT  ACCOUNTED  FOR  PREVIOUSLY 
00  4  I*1»NUMAIM 

. . DO  4  J*1 »  NGPW  ’  .  . 

X*J 

SECFRACI,  J)  =C  CX+FRAI)  *STEP*TOWIO*OAMRAO-OAHLOC Cl»  l/STADEV 
IFCCCX+FRAI) *STEP*CAMRAO*TOWIOl •GT.ERMIO) 

CSECFRA  C I, J) *  CERUNWI-OAFLOCCI J  )  /S  TAOEV 
4  CONTINUE 

00  6  K*1,NUNAIH  ‘  '  . .  . .  ~ 

DO  6  1*1, NGPW 

.  J*NS3G*I-1  . .  . ‘  .  . . . 

CALL  HONCRC  SECBON CK,J)  ,  BO  NOD 

CALL  MONORC  SECFRA  C K, I)  »BONOR>  “ . . 

PROF AC CK  ,1) *BONOR-BONO l 

6  —  CONTINUE  .  *  . 

Z  FORH  AN  30  ARRAY  OF  THESE  PROBABILITIES  WITH  1D*GR0UP  NUMBER 

C  2D*AI H  POINT  NUMBER,  3D=SECTI0N  NUMBER  CPROSPGI 

00  3  J*  1 ,  NGPW 
II=NSPG+ J-l 
00  3  L*l»  NUMAIM 
00  3  I*J,II 
K*I-J^1 

IFCK.LT.NSPGI  PROSPGC  J,L,  K) =  PROSEC(L» II 
IF  CK»  EQ»  NS°G)  °ROSPG  <J,L,K)*PROFACCL»J) 

3  CONTINUE 

RETURN 
ENO 
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Subroutine  UNION 


SUBROUTINE  UNION! NUM80M,N UMAIM, EPSILN, TOTAL) 


THIS  SUBROUTINE  (UNION)  USES  THE  EVENT  COMPOSITION 
METHOD  TO  CALCULATE  THE  PROBABILITY  OF  NOT  COMPLETING 
THE  CUT.  ONE  MINUS  THIS  NUMBER  (TOTAL!  THEN  GIVES 
THE  PROBABILITY  OF  MAKING  THE  CUT  WITH  THE  GIVEN  AIM  POINT 
STRATEGY. 

THE  FOLLOWING  VARIABLES  ARE* 

1.  T=AN  ARRAY  IN  WHICH  THE  SECTIONS  IN  A  GROUP  ARE  ADDEO 
U®  TO  DETERMINE  THE  PROBABILITY  FOR  THAT  GROUP 

2.  IMDX  =  AN  ARRAY  WHICH  IS  USEO  AS  A  POINTER  TO  DETERMINE 
WHICH  GROUP  IS  BEING  USEO  AND  ALSO  HOW  MUCH  OF  THAT 
GROUP  IS  TO  BE  USEO  IN  THE  CALCULATION 

3.  SUM* AN  ARRAY  WHICH  IS  COMPRISED  OF  THE  PROBABILITIES 
OF  THE  COMBINATIONS  OF  THE  GROUPING*  I.E.  SUM(l)  IS 
THE  SUM  OF  THE  PROBABILITIES  OF  ALL  GROUPING  TAKEN 

'*  ■'  ONE  AT  A  TIME.  SUM  (21  IS  THE  SUM  OF  THE  INTERSECTIONS 
OF  ALL  GROUPINGS  TAKEN  TWO  AT  A  TIME.  ETC. 

4.  EPSILN=A  CONVERGENCE  CRITERIA  ON  THE  CHANGES  IN  THE 
PROBABILITY  OF  MAKING  THE  CUT. 

5.  T3TSUM=THE  SUMMATION  OF  ALL  SUMCI! 


C 


DIMENSION  NUMBOM( 20) , INOX (?0)«T(20)«  SUM! 2  0) 
COMMON  PROSPGC 12, 20*201 
COMMON/ 1  NOICE/ INOX 

COMMON/FIXEO/ RUN  WI C*  0 AMRAO.S  TAOE V.TO WIO*  PROCUT 
COMMON/GRP/ NSPG.NSPW.NGPW,  STEP, FRAI,  GRUNWI,ERWIO 


C 
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TOTSUM  =  0 • 

00  101  K=l»  NGPH 
SUH(K) =0. 

K2=K*1 

00  11  11  =  1*  K2 
INOXC 1 1) =11 
CONTINUE 

INOXC  K )  =  INDX  <  K)  -1 
KK=K 

INOXC  KKi  =INOXCKK>  *1 
J=NGPM-( K-KKI 

IFCINOXC KK) .GT.J)  GO  TO  100 
IFCKK.EQ.K)  GO  TO  20 
M*  KK*  1 

00  15  KKK=M*  K 

INOXCKKK1=INOXCKK> ♦CKKK-KK) 

CONTINUE 

KK=K 

PR  00= 1  • 

00  30  L= 1* NUN AIM 
T(K)  =  0. 

00  25  I=1,NSPG 

TIKI  =  TCK)  *PROSPGC  INOXC  HO  ,L.I) 

CONTINUE 

s*i-rc  ki 

NM=K-1 

IFCMM.EQ.OJ  GO  TO  29 
00  2F  IK=1*  MM 
H1=IK*1 

IN=IN0XCM1) -INOXCIK) 

LIH=HINO  CNSPG*  IN) 

TCIK)=0. 

00  28  1=1, LIM 

T  C  IXI  *T  C  IK)  ♦PROSPG  (INOXCIK),  L»  I) 

CONTINUE 

S*S-TCIK) 

CONTINUE 

CONTINUE 

PROO*  PROC*C  S**NUMBOHCL )l 
CONTINUE 

SUHCK)  *SUHC  K)  ♦PROO 
GO  TO  10 


200 

101 

300 


l^i^s  Klf  *  1 

IF(KK.GT.O)  GO  TO  10 
B*KM 
2*8/2 

J5*<K*ll/2 
A*Z-J5 
J6=0 
J7=J6-  1 

IFtA.GT.  nSUH(K)*  J7*SUHCK> 
TOTSJM*TOTSUM  ♦SUM ( K I 
IF(K.EQ.l)  GO  TO  200 
U*ABS(SUKCK-m 
Vs ABS ( SUM (K ) ) 

H=U-V 

TOLAN*  A8  S(M ) 

IF(TOLAM.LT.EPSILN)  GO  TO  300 

CONTINUE 

CONTINUE 

TOTAL=TOTSUM 

RETURN 

ENO 
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Subroutine  SIM 


SUBROUTINE  SIM(NUMAIM*  AINLOC  t NUMB OM* TOTAL) 


C 


THIS  SUBROUTINE  ISIM)  USES  A  MONTE  CARLO  SIMULATION 
TO  DETERMINE  THE  PROBABILITY  OF  NOT  MAKING  THE  CUT. 

THIS  SUBROUTINE  IS  USED  WHENEVER  THE  NGPW  EXCEEDS  12. 

IF  NGPW  EXCEEDS  12  THIS  SUBROUTINE  IS  -QUICKER-  THAN 
UNION. 

THE  FOLLOWING  VARIABLES  ARE « 

1.  RANF(DUM)=A  RANDUM  NUMBER  GENERATOR 

2.  G3NML*IMSL  SUBROUTINE  WHICH  RETURNS  A  VECTOR  OF 
NORMAL  (0,1)  STANDARD  DEVIATES. 

3.  VSTRAsIMSL  SUBROUTINE  WHICH  SORTS  A  VECTOR  AND 
PLACES  THE  VALUES  IN  ASCENDING  ORDER.  THE  FIRST 
NUMBER  IN  THE  RETURNED  VECTOR  IS  THE  SMALLEST. 

4.  NIT  ERA=THE  NUMBER  OF  ITERATIONS  FOR  THIS  SIMULATION 


DIMENSION  NUMBOM  ( 20) , A  JML  OCC  20 ) 

DIMENSION  HIT (501 *R1( 1  ) 

COMMON /SIMS /H IT* R1 

COMMON/F IXED/RUNWID, OAWRA O.STA DEV ,TOHIO, PROCUT 
COMMON/ ITER AT/NITERA 
DOUBLE  PRECISION  OSEEO 


•  00  100  I»l»  50 
HIT<I)=0.0 
100  CONTINUE 

E*1000  00  000  • 

M*0 

C*0. 

C*C-3AMRAD 
0*0AMRAD*RU  NWID 
CLEAR=2.*OAMR«D*TOWIO 
ECLEAR=OAMRAO*TOWIO 
NUMBs 0 

00  1  1*1 »NUMA IN 
NUMB* NUMB+NUM BOM (  I) 

1  CONTINUE 
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1=0 

N=NITE  RA 

2  CONTINUE 
1=1*1 

IF(I.EQ.N*1)  CO  TO  30 
X=RANF  (OUM) 

IX=X*E 

IF ( IX •  EQ.  0)  IX=585285124 

OSEEO=IX 

N2=0 

DO  3  J  =  1  ,NUMAIM 
Ni=NUMBOM(J> 

00  3  J1=1,N1 
N2=N2*I 
12=1 

CALL  GGNML( DSEEO, 12, Rl) 

HIT(N2)  =AIMLOC( J>  *STADEV*R1I 12) 

3  CONTINUE 

CALL  VSRTAC HIT, NUMB) 

OOES  A  WEAPON  LAND  ON  THE  EFFECTIVE  RUNWAY 


N4=0 

N5=0 

N6=0 

DO  8  L*1.NUMB 
IF<HIT(L).LT.C)  N  4=N4*  1 
IFCHIT(L).GT.O)  N5=N5*1 
8  '  CONTINUE 

IFIN4. EQ.NUMB)  GO  TO  20  _ 

IFIN5. EQ.NUMB)  GO  TO  20 
N6=N4*N5 

’* IFCN6. EQ.NUMB)  GO  TO  20 
N4=N4*1  _ 

C  IS  THE  FIRST  WEAPON  CLOSE  TO  THE  RUNWAY  EOGE 


OIS=0. 

. . OIS=0.-H  IT<N4) 

IF(ABS (OIS) .GT.ECLEAR)  GO  TO  20 

C  IS  THE  "LAST"  WEAPON  CLOSE  TO  THE  RUNWAY  EOGE 

OIS*0. 

'  N9=N'JM8-N5 
0IS=RUNWI0-HIT(N9) 

IF(A3SCDIS)  .GT.ECLEAR)  GO  TO  20 


C  CHECK  ON  THE  OISTANCE  BETWEEN  ADJACENT  IMPACT  POINTS  ON  THE  EFFECTIVE 
C  RUNWAY 

DIS=0. 

IF ( (N4  *N5) • EQ. (NUMB-1 ) )  GO  TO  4 

N8=N JM  B-N5- 1 

IF(N4.EQ.N8)  GO  TO  4 

IF(N3 . EQ  «  0)  0IS=HIT<N4)-HITIN9) 

IF(N3.EQ.O)  GO  TO  4 
00  1  J=N4,N8 
DISMITC  JHI-HITC  J) 

IF(A9S  (OISI  .GT. CLEAR)  GO  TO  20 
7  CONTINUE 

4  CONTINUE 

IF(ABS(OIS) .GT. CLEAR)  GO  TO  20 
GO  TO  2 
20  CONTINUE 
Z=0. 

M=M*1 
GO  TO  2 
30  CONTINUE 
A=M 

total=a/b 

RETURN 

ENO 
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Subroutine  RUNPRO 


SUBROUTINE  RUNPROI Y, NCUT,  BNOLOW, BNDUP» 


C 


THIS  SUBROUTINE  CRUNPRC)  CALCULATES  THE  INTERVAL  FOR  THE  PROBABILITY 
OF  DESTROYING  THE  RUNWAT  AT  A  CONFIDENCE  LEVEL  OF  .99 
THE  FOLLOWING  VARIABLES  ARE* 

1.  E  S  T  X  =T  H  £  ESTIMATE  OF  THE  PROBABILITY  OF  DESTROYING  THE  RUNWAY 

2.  EVA  FY=THE  ESTIMATE  OF  THE  VARIANCE  IN  HEAN  FOR  ONE  CUT 

3.  ESTYSQ=THE  MEAN  SGUAREO  OF  A  CUT 

k,  EVARX=  ESTIMATE  OF  THE  VARIANCE  FOR  THE  RUNWAY  DESTRUCTION 

5.  STANDX=THE  STANDARD  DEVIATION  OF  THE  RUNWAY  DESTRUCTION 

6.  BNOLOW  =  THE  LOWER  EOUNO  ON  THE  INTERVAL 

7.  BNOUP=THE  UPPER  BOUND  ON  THE  INTERVAL 


C 


COMMOM/ITERAT/NrTERA 


C 


XlsNITERA 
'  X2*NCUT 
ESTX=V**X2 
EVARY*Y*  C l.-YI  /XI 
ESTYSQ=Y**2. 

EVARX= ( ( EVARY+ESTYSQ) **X2) -< ESTX **2. > 

STANDX*EVARX**.5 

BNOLOW*ESTX-3.*STANOX 

BNOUP*ESTX+3 .*STANOX 

RETURN 

END 
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APPENDIX  D 

Derivation  of  the  Mean  and  Variance 
of  P,  the  Estimated  Probability  of 
Destroying  the  Runway  With  the 


"Best"  Aim  Point  Strate 


r 


Derivation  of  the  Mean  and  Variance 

A 

of  P,  the  Estimated  Probability  of 
Destroying-  the  Runway  With  the 
"Best"  Aim  Point  Strategy 

One  of  the  basic  assumptions  used  in  this  model  is 
that  the  overall  problem  of  destroying  the  runway  can  be  ex¬ 
pressed  as  N  identical  and  independent  subproblems  (see  Note 
on  page  8) .  One  also  knows  the  distribution  of  the  estimates 
for  the  actual  value  of  p,  the  probability  of  making  a  cut 
with  a  given  aim  point  strategy.  If  one  lets:  1)  Y  be  the 
random  variable  which  denotes  the  estimates  of  p;  2)  Fy(y) 
is  the  distribution  of  Y;  and  3)  Y^,  Y2,  .  ..,  YN  are  i.i.d. 
from  Fy(y) .  Thus, 

E[Y!Y2. . .Yn]  =  E[Y1]E[Y2]...E[YN]  CD-I) 

where  E[Y]  is  the  expectation  of  the  random  variable  Y  (Ref  11 
179)  . 

Since  the  subproblems  are  identical,  i.e.  E[Y^]  = 

E[Y2]  *  . . .  *  E[Yn],  Eq  (D-l)  becomes 

E[Y1Y2--*YN1  =  CH[Y1])N  (D-2) 

Also,  according  to  Theorem  3.4  (Ref  11:  90) 

V[A]  -  E[A2]  -  y2  (D-3) 
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where  V[A]  is  the  variance  of  the  random  variable  A  and  y  is 
the  E  [A] .  Thus , 

V[Y1Y2'-*YN]  =  E[(V?...Yn)2]  ‘  (E(Y1Y2...Yn])2  (D-4) 

Coupling  the  assumption  that  the  Ys  are  independent,  identi¬ 
cally  distributed  random  variables  with  Eq  (D-3) 

V[Y1Y2...Yn]  =  (E[(Y1)2])N  -  (E  [Yx ] ) 2N  (D-5) 

According  to  Mendenhall  and  Scheaffer,  Y^  can  be  viewed 
as  asymptotically  normal  with  mean  p  and  variance  p(l-p)/M, 
where  M  is  the  number  of  simulated  attacks  (Ref  11:  256). 

Thus , 

E [ Y  x ]  =  p  (D-6) 

where  p  is  the  point  estimate  of  PC  for  the  "best"  aim  point 
strategy,  and 

V[YJ  =  (D.7) 

By  using  Eq  (D-4) 

E[(YX)2]  =  +  (p)2  (D-8) 

Thus 

vtYiY2...YN]  =  (21^-^-+  (P)2)N  -  (P)2N  (D-9) 

To  construct  a  99  percent  confidence  level  on  an  in- 

/v 

terval  estimate  of  P,  one  can  use  ±3o,  the  standard  deviation, 

A 

about  the  point  estimate  P.  To  accomplish  this  one  needs, 

£  “  (P)N  (D- 10) 
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and 


A 

O 


C(E£iiEl.  -  *2N)1/2 


(D- 11) 


(See  Appendix  C,  Subroutine  RUN?RO  for  the  coding  of  this 
derivation. ) 


i 
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A  generalized  computer  model  lor  determining  where  to 
target  conventional  weapons  against  a  non- reinforced  concrete 
runway  is  developed.  This  development  consists  of  taking  the 
overall  problem,  runway  destruction,  and  dividing  it  into  a 
number  of  identical,  stochastically  independent  subprob lems . 
Then  one  of  these  subproblems  is  used  to  determine  an  aim 
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tpoint  strategy  which  reduces  the  number  of  weapons  required 
to  achieve  the  level  of  runway  destruction  desired. 

In  order  to  solve  the  subproblem,  a  search  algorithm 
is  developed.  Along  with  this  search  algorithm,  two  methods 
for  approximatin'.;  the  level  of  destruction,  given  an  aim 
point  strategy,  are  investigated.  Finally,  the  entire  model 
is  verified  and  "validated"  hy  a  number  of  statistical  tests 
Also  a  comparison  of  this  model  to  the  Airfield  Assessment 
Program,  an  Air  lorce  Armament  laboratory  model  ,  is  provided 
Using  the  model's  results  lor  the  data  investigated,  an 
equation  relating  circular  error  probable  and  minimum  launch 
window's  width  to  the  number  of  weapons  employed  is  derived. 
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